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Chapter  I 
Introduction 


effort  made  for  this  study  contract  is  grouped  into 
two  major  categories: 

(1)  Experimental  and  theoretical  study  of  the  infrared  emission 
generated  in  a  low  pressure  gaseous  glow  discharge,  and 

(2)  Processing  of  SCRIBE  [Stratospheric  Cryogenic 
Interferometer  Balloon  Experiment]  data  . 

fx 


The  first  study  was  focused  on  the  laboratory  study  of 
the  infrared  emission  observable  in  the  atmosphere  above  100  km, 
while  the  second  was  on  the  field  data  of  the  infrared  emission 
which  was  observed  in  the  lower  atmosphere  below  30  km  using  a 
balloon-borne  cryogenic  interferometer  spectrometer.  Even  though 
these  studies  are  concerned  with  the  infrared  emission  of  the 
atmosphere,  their  characteristics  are  quite  different  .  It  would 
best  be  discussed  separately.  In  this  report,  the  study  of  the 
low-pressure  glow  discharge  source  will  be  discussed  in  Chapter 
II.  The  discussion  of  SCRIBE  data  will  be  made  in  Chapter  III. 

During  the  contract  period,  we  produced  twelve  papers, 
mostly  published  in  open  literature.  Our  Scientific  Report  I  was 
produced  by  assembling  three  papers  out  of  these  twelve  works. 
This  report  is  accompanied  by  most.  of  our  publications,  eight 


papers. 
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Even  though  Sakai  authored  this  report  by  taking  a  chief 
responsibility  of  assembling  the  works  accomplished  under  this 
contract,  his  contributions  to  those  published  works  are  limited. 

His  name,  which  appears  on  the  front  page  of  this  report,  should 
be  considered  as  to  signify  his  editor's  role  for  the  report. 

Many  people  made  contributions  leading  to  joint  publication. 

Some  of  them  were  either  principally  or  partially  supported  under 
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instrumental  in  completing  those  published  works. 
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Mark  P.  Esplin  and  Hajime  Sakai,  "High  Temperature  Absorption 
Spectrum  of  4.3  Micron  C02,"  presented  at  OSA  Topical  Meeting 
on  Spectroscopy  in  Support  of  Atmospheric  Measurements, 

Sarasota  FL,  November  1980. 

H.  Sakai,  "The  w^^-a^iT^  (0-0)  Band  of  N2,"  presented  at  OSU 
36th  Annual  Symposium  on  Molecular  Spectroscopy,  June  1981. 

E.S.  Chang,  "Theory  of  IR  Spectra  from  Rydberg-  Rydberg 
Non-penetrating  State  Transitions,"  and 

L.S.  Mayants,  "New  Rotation-Vibration  Hamiltonian,"  presented 
at  OSU  Molecular  Spectroscopy  Symposium,  June  1982. 

H.  Sakai  and  G.  Vanasse,  "Stratospheric  Infrared  Emission 
Observed  by  Balloon-Borne  Cryogenic  Interferometer 
Spectrometer,"  OSA  Annual  Meeeting,  October  1982. 

H.  Sakai,  "Molecular  Constants  of  12CL602  Bands  in  1900  cm  1  ^ 
2150  cm  1 ,"  and 

E.S.  Chang,  "Identification  of  the  nf  Complexes  in  N 2  by  UV 
Absorption";  and 

E.S.  Chang,  "Absorption  of  the  5g-4f  Emission  Lines  in  Hg," 
presented  at  OSU  Molecular  Spectroscopy  Symposium,  June  1983* 


Chapter  II 


Study  of  the  infrared  emission  generated 
in  a  low-pressure  gaseous  glow  discharge 

Our  study  was  conducted  using  a  12-meter-long  glow 
discharge  column  shown  in  Figure  1.  In  the  early  part  of  our 
study  ,  we  found  that  the  source  was  good  in  generating  various 
species  involved  in  the  infrared  emission  of  the  upper 

atmosphere.  However,  their  spectral  characteristics  were  rather 
uncontrollable;  we  were  not  able  to  control  their  temporal 
behavior  concerning  the  excitation  or  the  relaxation.  Their 

observed  radiative  temperatures  were  so  extraordinary  that  the 
spectral  data  provided  no  clue  for  analysing  what  was  happening 
in  the  glow  discharge  process.  We  decided  to  develop  the 

time-resolved  technique  for  making  a  detailed  observation  of  the 
infrared  emission  at  various  phases  of  excitation-relaxation 

process  in  our  glow  discharge  source. 

We  set  the  primary  experimental  objective  to  detect  the 
atmospheric  infrared  emission  species  as  many  as  possible,  thus 
allowing  the  chemistry  processes  in  the  source  to  occur  in  a 
maximum  extent.  The  flow  rate  of  molecules  in  the  discharge 

region  was  very  3low.  A  molecule  injected  in  the  source  was 
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subject  to  many  discharge  cycles,  in  an  order  of  10  ,  before  it 


was  pumped  out  from  the  source  chamber.  The  discharge  cycle  was 
set  to  be  as  slow  as  possible  for  achieving  a  stable 
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excitation-relaxation  cycle  in  the  emission.  As  a  result,  we 
were  able  to  pick  up  various  species,  and  at  the  same  time  we 
established  no  control  over  their  excitation-relaxation  evolution 
sequence . 

The  glow-discharge  source  was  activated  by  applying  an 
ac  house-current  of  60  Hz.  The  voltage  at  the  power  supply  was 
varied  up  to  1000  V  r.m.s,  or  1400  V  peak.  The  ballast 
resistance,  connected  in  series  with  the  discharge  electrode, 
somewhat  regulated  the  discharge-voltage  to  a  constant  value  of 
600  V.  This  figure  was  observed  to  be  held  constant,  independent 
of  the  supply  voltage.  The  current  was  varied  by  setting-  the 
supply-voltage  to  a  desired  value.  Thus  by  setting  the 
supply-voltage  to  a  higher  value,  we  were  able  to  feed  a  heavier 
discharge  current.  There  was  a  distinctively  different  current 
behavior  between  the  two  opposite  polarity  phases;  i.e.,  between 
the  phase  when  the  central  electrode  was  positively  charged  and 
that  negatively  charged.  The  single-polarity  operation  led  to 
unstable  discharge,  unfit  to  our  experimental  measurement.  The 
central  electrode  had  to  be  alternately  charged  with  opposite 
polarity  in  order  to  achieve  a  stable  glow  discharge  in  the 
chamber . 

The  technique  of  time-resolved  Fourier  spectroscopy  is 
outlined  in  one  of  our  papers  [  Barowy  and  Sakai  (1984)  ].  One 
of  the  experimental  requirements  imposed  on  us  was  to  gain  a 


reasonable  signal-to-noise  ratio  in  the  observation  of  a 
relatively  weak  emission  source.  Since  our  source  emission  was 
generated  in  a  repetitive  pattern,  our  experiment  was  best  done 
by  implementing  a  step-and-hold  scanning  to  the  interferometer 
drive.  That  is  to  say,  the  interferometer  was  held  at  a  given 
optical-path-difference  position  for  an  indefinite  duration  to 
accumulate  a  desired  amount  of  the  signal  until  a  reasonable 
signal-to-noise  ratio  was  achieved.  Our  scheme  was  more 
versatile  than  those  developed  by  the  French  group,  when  the 
system  was  applied  to  the  time-resolved  measurement. 

For  maintaining  a  technical  simplicity  to  our 
interferometer  instrumentation,  we  decided  to  use  the  intensity 
of  the  laser-fringe  signal  as  the  servo-control  signal. The 
alignment  stability  in  the  interferometer  circuit,  which  was  a 
paramount  importance  for  such  a  servo-control  scheme,  was 
provided  by  an  automatic  optical  retroreflecting 
tilt-compensating  scheme.  We  adopted  a  pair  of  the  readily 
available  lightweight  hollow  corner-reflectors,  replacing  a  set 
of  the  plain  mirrors  in  both  arras  of  the  interferometer.  At  the 
same  time  a  silicone  fluid  of  higher  viscosity  than  the  one 
supplied  in  the  original  interferometer  was  tried,  hoping  to 
attain  a  quicker  dissipation  of  the  kinetic  energy  of  the  movable 
carriage.  The  interferometer  drive  mechanism  designed  with  a 
large  damping  term  is  to  achieve  a  fast-acting  stepping  drive 
motion.  At  present  ,  we  find  the  new  damping  fluid  and  the 
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existing  dash  pot  structure  in  the  original  design  provides  an 
adequate  stepping  speed  of  approximately  10  steps  per  second. 
The  present  operation  does  not  require  the  stepping  speed  of  this 
magnitude,  as  the  measurements  deal  with  a  weak  source.  The 
drive  motor  was  energized  by  a  series  of  the  fast  impulses.  The 
control  electronics  and  other  pertinent  servo-logics  are 
described  in  our  paper  cited  above. 

Replacement  of  the  plain  mirrors  with  the  hollow 
retro-reflecting  corner  mirrors  in  the  interferometer  arms 
lowered  the  mechanical  precision  required  on  the  way  system  of 
the  interferometer  drive.  The  movable  carriage  wa3  held  much 
more  loosely  than  the  conventional  design,  resulting  in  a 
substantial  reduction  of  the  harmful  static  friction  in  its 
motion. 

The  data  shown  in  Figure  2  are  the  time-resolved 
spectra  of  the  glow-discharge  emission.  The  data  shown  in 

Figure  3  are  those  of  the  0 2  emission.  It  is  noticeable  that 
the  emission  generated  in  an  electric  polarity  phase  differs 
vastly  from  that  in  another  phase.  The  electric  field  at  the 
central  electrode  of  our  source  configuration  is  approximately 
13.5  times  stronger  than  that  at  the  outer  electrode.  The 
electronic  emission  by  the  field  is  likely  to  occur  only  when  the 
central  electrode  is  negatively  charged  [the  negative  phase]. 
The  electronic  emission  from  the  outer  electrode  is  less  likely 
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even  when  that  electrode  is  negatively  charged  [the  positive 
phase].  The  current  observed  during  the  positive  phase  is  very 
likely  carried  by  the  positive  ions,  which  have  the  mobility 
approximately  100  times  slower  than  the  electrons.  As  a 
consequence,  the  electric  conduction  during  the  negative  phase  is 
done  by  both  the  electrons  and  the  positive  ions,  while  the 
electrons  are  depleted  from  the  conduction  during  the  positive 
phase.  Meanwhile,  the  current  magnitudes  are  observed  to  be 
equal  for  both  polarity  phases.  This  does  not  mean  that  the 
current  carriers  in  the  discharge  chamber  are  equalized  for  both 
phases.  It  indicates  that  the  total  amount  of  the  current 
leaving  the  discharge  chamber  remains  the  same  for  both  phases. 
A  majority  of  the  carriers  moves  toward  the  outer  electrode 
during  the  positive  phase,  while  during  another  phase  two  kinds 
of  the  carriers  move  in  .  opposite  directions.  Those  species 
excited  by  the  electron  collision  are  very  likely  observable  only 
during  the  negative  phase.  The  time-resolved  measurement  is  very 
suitable  to  observe  the  difference  in  the  radiative  emission 
between  these  two  phases.  It  provides  a  very  powerful  means  for 
revealing  an  intricate  change  in  the  spectral  characteristics  at 
various  stages  of  the  excitation-relaxation  sequence.  The 
time-resolved  data  shown  in  Figures  2  and  3  demonstrate  the 
capability  provided  by  the  technique  in  following  the  temporal 
sequence  of  each  spectral  system.  For  our  primary  goal  of 
detecting  a  particular  atmospheric  species,  we  can  find  the  most 
favorable  time  phase  for  the  observation.  The  spectrum  of  NO 


shown  in  Figure  1*  gives  another  example  of  the  power  of  the 
time-resolved  Fourier  spectroscopy  technique. 

During  the  course  of  our  spectral  analysis  effort,  we 
came  to  realize  the  structure  observed  at  a  vicinity  of  2500  cm 
in  various  spectral  data  can  be  interpreted  as  the  4f-5g 
transition  of  either  atomic  or  molecular  species.  We  had  a  good 
success  in  identifying  the  4f-5g  transition  lines  of  the  atomic 
nitrogen  for  both  the  line  positions  and  the  line  intensity. 
Subsequently  our  theoretical  effort  was  extended  principally  by 
E.  Chang  and  several  co-operative  efforts  brought  a  fruitful 


development. 
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Chapter  III 


SCRIBE  Data  Processing 

The  work  performed  under  the  present  contract  for  the 
AFGL  SCRIBE  program  was  an  ongoing  effort.  The  program  itself 
has  been  described  in  several  of  our  publications  as  well  as  in 
our  previous  reports.  Our  participation  may  be  summarized  into 
three  major  parts: 

(1)  Extraction  of  the  interferogram  data  from  the  PCM  telemetry 
data  tape; 

(2)  Processing  of  the  interferogram  data  for  the  spectral 
retrieval;  and 

(3)  Analysis  of  the  obtained  atmospheric  emission  data. 

Our  Scientific  Report  II  of  this  contract  was  assembled  on  the 
SCRIBE  II  flight  data  [  Oct.  1982  3.  The  parts  (2)  and  (3)  have 
been  documented  in  our  previous  reports.  The  following 
description  is  on  the  part  (1),  the  extraction  of  the 
interferogram  data  from  the  PCM  telemetry  data  tape. 

Introduction 

The  stratospheric  cryogenic  interferometer  balloon 
program  is  a  balloon-borne  experiment,  intending  to  obtain  the 
stratospheric  emission  data  observed  along  the  horizontal  line  of 
sight.  Four  groups  are  involved  from  the  onset  of  the  program, 
AFGL,  Idealab  Inc.,  University  of  Denver  and  University  of 


Massachusetts. 
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The  obtained  data  during  the  balloon  flight  are 
transmitted  to  the  ground  station  via  a  radio  communication  link. 
The  instrumentation  for  this  radio  link  was  developed  under  a 
direction  of  the  AFGL  Balloon  Branch  which  is  in  charge  of  the 
balloon  flight  operation. 

PCM  Telemetry  Signal 

The  basic  form  of  the  telemetry  signal  is  a  72-bit 
sequence  in  a  binary  bipolar  non— return— zero— level  [NRZL]  format. 
The  signal  is  constructed  in  a  form  of  the  pulse-coded  modulation 
signal  [PCM].  All  the  pertinent  data  for  the  flight-experiment 
including  the  digitized  interferogram  data  are  processed  by  the 
on-board  PCM  encoder  into  a  72-bit  frame  pattern.  Each  frame  is 
made  of  a  20-bit  synchronization  identification  word,  a  4-bit 
subframe  identification  counter  number,  a  16-bit  interferogram 
data,  and  four  8-bit  housekeeping  data.  The  synchronization  ID 
code  shown  in  Table  I  provides  a  key  for  identifying  the  basic 
PCM  clock  rate  and  the  frame  length.  Ihe  interferogram  data  are 
a  12-bit  output  of  the  detector  amplifier  signal  and  a  4-bit  word 
which  indicates  the  status  regarding  its  validity,  the 
interferometer  drive  status,  and  the  amplifier  gain  status. 

The  telemetry  data  are  transmitted  from  the  flight 
package  at  a  rate  of  800  K  bit  per  second,  and  they  are  recorded 
on  the  analog  magnetic  tapes  at  the  ground  station.  The  spectral 
recovery  process  in  Fourier  spectroscopy  requires  all  the 
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inter ferogram  data  taken  at  every  sampling  position.  The 
extraction  process  of  the  inter ferogram  data  ,  therefore,  allows 
no  failure  in  recovering  them  from  the  recorded  telemetry  data 
which  are  subjected  to  various  electrical  interferences.  A 
considerable  time  and  effort  were  spent  by  us  for  achieving  an 
improved  reliability  of  recovering  the  interferogram  data  from 
the  telemetry  record. 

General  Layout  of  Processing  Scheme 

The  data  extraction  scheme  currently  used  at  our 
facility  is  shown  in  Figure  5  .  The  PCM  telemetry  record  on  the 
magnetic  tape  was  played  back  by  using  the  Honeywell  101  tape 
recorder.  The  recovered  PCM  signal  was  sent  to  the  PCM 
decommutator,  which  produced  two  outputs,  a  parallel  16-bit 
signal  and  a  clock  signal  synchronous  to  the  incoming  PCM  signal. 
The  home-built  interface  electronics  between  the  PCM  decommutator 
and  the  PDP  11/20  computer  provided  a  temporary  storage  of  a 
proper  16-bit  word  extracted  from  the  serial  72-bit  string  of  the 
telemetry  signal,  and  it  sent  a  16-bit  output  to  the  DR11-B  I/O 
interface  of  the  PDP  computer.  The  data  transferred  into  the 
central  memory  of  the  PDP  were  eventually  stored  on  the  disk 
memory  [RK1],  and  at  the  end  they  were  copied  onto  the  magnetic 
tape  using  the  PIP  [Peripheral  Interchange  Program]  of  the  RT11 
software  package. 
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Block  Diagram  of  PCM  Data  Extraction  Scheme 


The  data  extraction  was  essentially  controlled  by  two 
factors,  a  signal-to-noise  ratio  of  the  played-back  PCM  signal, 
and  a  timing  jitter  contained  in  it.  The  computer  simulation 
test  showed  that  with  the  specified  bandwidth  in  the  PCM 
electronics  train  the  timing  jitter  would  not  exceed  a  figure 
more  than  a  128-th  of  the  basic  PCM  pulse  width,  even  if  the 
harmonic  distortion  in  the  played-back  signal  appeared  to  be 
rather  bad.  We  found  that  a  pulse  shaping  network  made  by  using 
the  broadband  operational  amplifier,  Tektronix  AM501,  worked 
adequately  for  our  setup.  The  PCM  decommutator  used  was  designed 
for  a  bit  rate  much  slower  than  our  800  K  bit  per  second,  and  it 
was  modified  in  a  patchwork  manner  until  we  obtained  the  latest 
data. By  taking  an  advantage  of  relatively  less  busy  time 
available  to  us  between  the  flights,  the  PCM  decommutator  was 
extensively  studied  and  the  modifications,  which  we  concluded  to 
be  necessary  from  the  study,  were  permanently  implemented. 

PCM  Decommutator 

The  original  design  of  our  PCM  decommutator  called  for  a 
rate  of  10  Kbit  per  second  ,  which  was  far  lower  than  our  800 
KBS. In  addition  its  synchronization  logic  assumed  an  absence  of  a 
long  non-changing  bit  pattern,  which  was  always  present  in  our 
PCM  signal.  We  concluded  that  an  extensive  redesign  of  some 
logics  was  necessary  to  achieve  improved  reliability  of  the 
operation. 
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The  basic  scheme  used  in  our  system  is  sketched  in 
Figure  6  .  The  transition  detector  provides  a  sharp  pulse  every 
time  when  the  PCM  signal  undergoes  an  amplitude  transition.  The 
SPD  phase  detector  operates  on  the  clock  signal  generated  at  the 
V CO  [voltage  controlled  oscillator]  timing  circuit  and  the  PCM 
transient  pulse;  every  time  when  these  pulses  arrive  within  a 
prescribed  time  interval,  the  SPD  triggers  the  storage  register 
to  pass  the  output  of  the  *J-bit  counter  to  the  digital-to-analog 
converter.  The  servo-error  signal  to  the  VCO  which  operates  at 
an  approximate  rate  of  16  times  the  basic  PCM  clock  rate,  is 
generated  only  at  an  instant  when  the  transition  detector 
provides  a  pulse  to  the  SPD.  A  lengthy  non-changing  bit  pattern 
which  is  present  in  our  PCM  signal  makes  the  servo-loop  in  the 
VCO  inactive,  and  thus  the  VCO  locks  its  clock  rate  at  the  same 
value.  Consequently  no  shifting  of  the  VCO  clock  rate  occurs  in 
such  a  pulse  train.  All  circuits  which  deal  with  fast  pulse 
train  are  constructed  using  the  fast  Schottky  series  IC.  No 
significant  time  delays  are  observed  in  our  circuit  when  we  apply 
the  circuit  for  our  800  KBS  PCM  signal. 
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Improved  Bit-Synchronization  Scheme 
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Abstract.  From  the  spectroscopic  data  of  non-penetrating  Rydberg  states  of  C  l,  N  t,  O  l 
and  Ne  l,  we  have  extracted  the  polarisability  and  the  quadruple  moment  of  the  ionic 
core.  The  values  so  obtained  from  different  Rydberg  nl  levels  are  mutually  consistent  to 
within  1%.  Comparison  with  the  values  from  ab  initio  calculation  shows  agreement  to  a 
few  per  cent.  Further,  some  empirical  relations  between  certain  Slater  integrals  are  shown 
to  be  the  direct  consequence  of  our  treatment. 


It  has  long  been  known  that  the  ionic  polarisability  of  close-shell  cores  can  be  extracted 
from  the  energy  levels  of  their  non-penetrating  Rydberg  states  (Edl6n  1964).  For  the 
case  of  the  two-electron  core,  Olme  (1969)  has  shown  that  the  values  of  the  polarisabil¬ 
ity  of  Li*,  Be2*  and  BJ*  so  obtained  agree  with  the  best  theoretical  values  to  better 
than  one  per  cent.  In  the  case  of  open-shell  cores,  Rydberg  states  are  usually  analysed 
in  the  ‘pair-coupling*  theory  (Eriksson  1956)  or  equivalently  in  the  yAT-coupling  scheme 
(Humphreys  et  al  1967).  Such  analyses  yield  the  average  energy  of  each  nl  level  (the 
centre  of  gravity),  the  core  fine-structure  splittings,  the  Slater  integral  F2(nt)  and  other 
parameters  which  are  usually  very  small.  Using  these  determined  quantities,  known 
states  can  be  theoretically  obtained  accurately  and  unknown  states  can  be  predicted. 
While  the  value  of  the  polarisability  had  been  determined  from  the  average  energy 
in  some  open-core  cases,  it  was  thought  that  its  interpretation  as  the  core  polarisability 
had  to  be  modified  (Litzen  1968).  Consequently,  no  attempt  was  made  to  compare 
its  value  with  the  ones  from  ab  initio  calculations. 

In  this  letter,  we  determine  the  polarisability  of  atomic  cores  with  the  configuration 
2 pq  from  the  average  energy,  and  show  that  it  agrees  with  theoretical  values  to  a  few 
per  cent.  More  importantly,  we  find  that  the  electric  quadrupole  moment  of  the  core 
can  be  obtained  directly  from  the  F2(nl)  integrals.  Finally,  the  two  empirical  relations, 
(i)  niF1(nl)  =  constant,  noted  by  Eriksson  and  Johansson  (1961),  and  (ii) 
F2(5g)/F2(4f)  =  0.0466,  noted  by  Eriksson  (1967),  are  derived  explicitly.  In  fact,  all 
F2(nl)  are  expressible  as  a  product  of  the  core's  quadrupole  moment  and  an  analytically 
known  proportionality  constant. 

The  system  under  investigation  consists  of  an  atomic  core  with  a  well  defined 
fine-structure  labelled  by  the  quantum  number  /  (/  =LC+Se)  and  a  non-penetrating 
Rydberg  electron  with  orbital  angular  momentum  /  (typically  /  ?  3).  The  interaction 
between  the  electron  and  the  core  is  not  represented  by  the  usual  quantum  defect 
(because  the  electron  does  not  penetrate  the  core),  but  by  long-range  interactions. 
The  two  dominant  interactions  are  well  known  to  be  the  polarisation  potential  and 
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the  electric  quadrupole  potential  given  in  au  by 
e2a  eO 

V  =  -^-~P2(cos6). 


(1) 


In  equation  (1),  a  is  the  isotropic  polarisability  and  Q  is  the  quadrupole  moment 
defined  by  (Fano  and  Racah  1959) 

0  =  -e(Lc\\Cm\\Lc)(r2).  (2) 

Treatment  of  equation  (1)  by  perturbation  theory,  given  elsewhere  (Chang  and 
Sakai  1982),  yields  results  mathematically  equivalent  to  pair-coupling  theory.  For 
non-penetrating  states,  the  Rydberg  electron  is  always  exterior  to  the  core  electrons. 
Therefore  in  the  Slater  integral,  r<  refers  only  to  the  core  (2p)  electrons  and  r>  to 
the  Rydberg  (nl)  electrons,  and  the  integral  is  then  given  by 

F\2x>,nl)  =  e\r\v{r\l.  (3) 

The  last  integral  is  given  by  Bethe  and  Salpeter  (1967), 

<r’3>"'V/(M)(/+ 1)  (4) 


where  (  is  the  net  charge  on  the  ionic  core.  With  the  definition  F2(2p,nl)  = 
F2( 2p,  nl)/D2(l)  where  D2( 3)  =  75  and  D2(4)  =  385,  it  is  evident  that 


n  *F2(nl)  = 


eVV)  2P 

ni+b(i+\)D2{iy 


(5) 


Clearly  the  right-hand  side  of  equation  (5)  is  a  constant  for  a  particular  atomic  species 
and  for  a  fixed  value  of  l,  as  required  by  the  first  empirical  observation  (Eriksson  and 
Johansson  1961).  Now  from  equations  (3)  and  (4),  the  ratio 


F2(2p,5g)/F2(2p,4f)  = 


43x3xjx4x£>2(3)  43 

53x4xfx5x£>2(4)  53xll 


0.046545 


which  is  the  second  empirical  rule  (Eriksson  1967). 

We  now  proceed  to  evaluate  the  core  properties,  a  and  <?,  from  spectroscopic 
data.  The  evaluation  of  a  follows  the  standard  procedure  (Edlen  1964).  The  energy 
shift  Ap  of  a  Rydberg  level  is  usually  already  found  or  can  be  obtained  from  the 
average  energy  E0  and  the  series  limit.  Then  Ap  =  ( *aP{n ,  /),  where 


#*(»,/) -<r  A. 


3n2-/(/  + 1) 

2  n’(l-bt(l  +  i)l(l  +  l) 


(6) 


is  tabulated  by  Edlen  (1964).  We  have  ignored  the  contribution  due  to  the  induced 
quadrupole  polarisability  because  a  non-adiabatic  correction  was  about  the  same  order 
of  magnitude  but  opposite  in  sign  (Deutsch  1970).  Table  1  shows  the  polarisabilities 
of  C+,  N*,  O*  and  Ne"  so  obtained  from  the  spectra  of  C  i,  N  i,  O  i  and  Ne  i.  Note 
that  for  each  ion  the  polarisability  extracted  from  different  Rydberg  levels  is  consistent 
to  within  1%.  Two  exceptions  are  easily  explained  by  the  original  spectral  lines  being 
blended,  and  therefore  inaccurate. 

As  indicated  in  equation  (2),  the  quadrupole  moment  contains  the  radial  matrix 
element  (r2)?p  which  in  turn  can  be  extracted  from  the  spectroscopic  constant  F2  as 
shown  in  equation  (3).  Table  2  shows  F2  and  (r2)2p  obtained  from  the  same  Rydberg 


Letter  to  the  Editor 


L651 


Table  1.  Polarisability  from  Rydberg  levels. 


Species 

Cl 

Ni 

Ol 

Nei 

nl  A p  (cm‘ 

‘)  a(flo) 

Ap  (cm* 

')  a(aj) 

Ap  (cm 

')  a(flo) 

Ap  (cm* 

')  a(aj) 

4f  23.22" 

5.69 

14.99" 

3.67 

10.42' 

2.55 

5.41h 

1.33 

5f  13.38" 

5.72 

8.59" 

3.67 

6.03' 

2.58 

3.11h 

1.33 

6f  8.20" 

5.72 

S.16d 

3.60 

3.7l' 

2.59 

1.90h 

1.33 

7f  5.34" 

5.73 

2.42' 

2.59 

1.23h 

1.32 

8f  3.81"b 

(5.97) 

0.84'’ 

1.32 

9f 

0.60h 

1.31 

5g 

2.05" 

3.68 

0.74h 

1.32 

6g 

0.97u 

(2.72) 

"Johansson  (196S). 

bSome  components  blended,  some  with  J  unidentified. 

'Eriksson  and  Johansson  (1961). 
dMcConkey  el  at  (1968). 

"Chang  and  Sakai  (1981). 

'Moore  (1976). 
‘Blended. 
hLitzen  (1968). 

Table  2.  (r2)np  from  Rydberg  levels. 

Species 

Cl 

Nl 

Ol 

Ne  I 

nl 

F2(cm  1 

)  (^2)2p(^o) 

F2  (cm-1) 

(f2)2p(flo) 

F2I  cm 

')  (r2>2p(ao) 

F2  (cm*1)  (r\„(a20) 

4f 

2.80" 

2.57 

1.965" 

1.805 

o' 

_ 

r71.42h/n5  1.0251 

5f 

1.45* 

2.60 

1.000" 

1.794 

o' 

_ 

6f 

0.84" 

2.60 

0.580d 

1.798 

o' 

— 

7f 

0.54* 

2.61 

o' 

— 

8f 

0.35" h 

(2.57) 

9f 

5g 

j  0.0948" 

1.871 

10.08667 

1.710 

6g 

o'.« 

— 

0.05 19h  1.024 

"  hSame  as  in  table  1 . 


levels  and  atomic  species  as  in  table  1.  Once  again,  the  consistency  is  within  1%.  In 
fact  for  Ne  i,  F2(nf)  fitted  the  n  1  relation  to  within  the  experimental  error  and  were 
not  given  separately  for  each  n  (Litzen  1 968).  In  the  case  of  N  I,  our  F2(5g)  is  obtained 
from  only  two  states  (the  unperturbed  K  =  2  and  6)  of  the  data  of  Chang  and  Sakai 
(1981).  This  value  differs  by  8%  from  the  theoretical  value  0.086  67  cm'1  of  Eriksson 
(1967).  A  more  reliable  value  will  be  available  when  our  data  is  fully  analysed  (Chang 
and  Sakai  1982).  Technically,  the  value  for  F2  or  (r:>2p  extracted  from  spectroscopic 
data  includes  higher-order  terms  than  just  the  first-order  energy.  Consequently,  the 
quadrupole  moment  so  obtained  is  the  effective  quadrupole  moment,  which  includes 
a  small  contribution  from  the  anisotropic  polarisability. 
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To  complete  the  determination  of  the  quadrupole  moment,  we  need  to  evaluate 
the  reduced  matrix  element  in  equation  (2).  For  the  2P  core,  Fano  and  Racah  (19S9) 
give 

("C'!’|,)-3(J  J  S’4 

The  other  configurations  of  interest  p2  and  p3  are  treated  in  Slater  (1960).  Obviously 
configurations  2p6~q  give  the  same  values  except  for  the  sign,  as  negative  electrons 
are  replaced  by  positive  holes.  These  results  are  displayed  in  table  3. 

Finally,  we  summarise  our  atomic  properties  by  showing  their  mean  values  with 
standard  deviations  and  compare  them  with  theoretical  values  in  table  4.  Our 
polarisabilities  are  in  excellent  agreement  with  the  polarised  pseudo-states  calculation 
of  Hibbert  et  al  (1977).  Similarly,  our  quadrupole  moments  agree  very  well  with 
theoretical  values,  even  though  those  of  Sen  (1979)  are  for  a  different  system  with 
an  additional  3s  electron. 

We  conclude  that  the  polarisability  and  the  quadrupole  moment  of  an  atomic  core 
can  be  extracted  from  non-penetrating  Rydberg  states  with  an  accuracy  of  1%.  These 
values  are  useful  in  the  calculation  of  photoionisation,  free-free  transitions  and  as 
tests  for  wavefunctions  and  computational  methods.  In  spectroscopy,  we  have  derived 
the  two  empirical  rules  observed  by  Eriksson  and  will  show  how  the  core  properties 
can  be  used  to  predict  other  Rydberg  states. 


P 


Table  3.  The  reduced  matrix  element  in  equation  (2)  -«(L||C121|lL). 


""  L 

Configuration"' — 

S 

P 

D 

np 

-'A 

"P! 

+eA 

np 

0 

np4 

-'A 

»P5 

+eVf 

Table  4.  Polarisability  and  quadrupole  moment. 

a(<*o) 

Q(eaJ) 

Present 

Theory 

Present 

Theory 

C‘ 

5.72  ±0.02 

5.62* 

-2.82  ±0.02 

-2.834, b  -2.804' 

N* 

3.66  ±0.04 

3.63* 

1.97  ±0.01 

2.024c 

O4 

2.58  ±0.02 

2.57‘ 

0 

0 

Ne4 

1.32±0.01 

1.27* 

1.122±0.001 

1.283c 

•Hibbert  et  al  (1977). 
"Shorer  (1982). 

•Sen  (1979)  actually  2p"  3*. 


We  thank  A  Dalgamo  for  useful  correspondence.  This  work  was  supported  by  Air 
Force  Contract  No  F19628-81-K-0007. 
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Theory  of  angular  distributions  of  electrons  resonantly  scattered  by  molecules. 
III.  The  broad  resonance  in  N2  and  Co  at  20  eV 

Edward  S.  Chang 

Department  of  Physics  and  Astronomy,  University  of  Massachusetts,  Amherst.  Massachusetts  01003 
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Theoretical  expressions  for  the  vibrational  excitation  differential  cross  section  (DCS)  are 
given  for  the  71f  and  the  2 2.  *  resonances.  They  fit  some  recent  data  on  N.  and  CO, 
respectively,  at  an  electron  impact  energy  of  about  20  eV.  These  fits  allow  us  to  obtain  indi¬ 
vidual  rotational-vibrational  DCS’s.  It  is  demonstrated  that  the  discrepancy  between  the 
DCS  data  from  different  experiments  is  due  to  varying  energy  resolution  resulting  in  selec¬ 
tive  partial  exclusion  of  some  rotational  branches.  Further,  theoretical  DCS's  for  each  rota¬ 
tional  branch  in  both  gases  are  found,  which  are  verifiable  with  the  new  experimental  tech¬ 
nique  used  to  separate  out  the  rotational  branches  at  2  eV. 


It  is  well  known  that  resonances  enhance  the  vi¬ 
brational  excitation  cross  sections  of  molecules  by 
electron  impact.1  In  such  instances,  the  angular  dis¬ 
tribution  of  the  scattering  electron  or  the  differential 
cross  section  (DCS)  can  often  be  predicted  without 
an  ab  initio  calculation.  Such  a  theory  was  given  in 
an  earlier  paper  referred  to  as  I,  and  an  application 
to  CO  at  1.8  eV  was  discussed1  in  II.  In  the  present 
work,  we  apply  the  theory  developed  in  I  to  the  vi¬ 
brational  and  rotational  excitation  of  N,  and  CO  at 
about  20  eV. 

The  enhancement  of  the  vibrational  cross  section 
in  Ni  was  first  observed  by  Pavlovic  et  al.,*  who 
placed  the  position  of  the  resonance  at  22  eV  with  a 
width  of  3  to  4  eV.  Calculations  based  on  the  multi¬ 
ple  scattering  model  by  Dehmer  et  a!.'  show  a  :X„k 
resonance  in  the  15-  to  35-eV  range.  They  recog¬ 
nized  that  the  identification  was  not  complete  until 
the  DCS  was  calculated  but  were  discouraged  by  the 
strong  dependence  of  the  observed  DCS  with  energy 
in  the  resonance  region,  refuting  the  postulate  of  a 
single  resonance.4  Fortunately,  more  recent  mea¬ 
surements  of  the  vibrational  DCS  in  CO  as  well  as 
in  N2  by  Tanaka6  7  and  his  collaborators,  and  in 
both  gases  by  Tronc  et  al.,*  showed  little  variation 
with  energy.  These  data  enable  us  to  fit  theoretical 
angular  distributions  obtained  by  an  extension  of  I. 
Thus  we  are  able  to  reject  other  possible  resonances 
such  as  ftr,  fb,  and  f6,  while  confirming  the  fa  res¬ 
onance  in  both  N2  and  CO.  However,  the  best  fit 
with  experiment  is  achieved  only  when  the  res¬ 
onance  in  N2  admits  some  p  component,  and 
resonance  in  CO  some  d  component. 

Further,  the  same  theory  yields  the  individual 
rotational-vibrational  DCS's  for  the  'If  and  the 
:X  +  resonance.  It  can  be  seen  that  the  theoretical 
DCS’s  for  the  higher  rotational  branches  are  very 
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different  from  the  (elastic)  Q  branch.  Therefore 
their  exclusion  through  higher  energy  resolution, 
e.g.,  30  meV,  will  result  in  serious  distortion  of  the 
measured  vibrational  DCS.  In  particular,  this  effect 
appears  to  account  for  the  very  different  DCS  of 
Pavlovic  et  al.*  from  the  other  two  experimental 
groups.  Independently,  the  theoretical  rotational- 
vibrational  DCS  can  be  directly  verified  with  the 
new  experimental  technique  of  Jung  el  al.?  as  their 
2-eV  data  has  confirmed  the  theory  of  the  2I1  reso¬ 
nance  discussed  in  II. 

In  I,  we  have  derived  an  expression  for  the  angu¬ 
lar  distribution  of  electrons  scattered  from  a  I-state 
molecule.  It  is  written  as  a  sum  over  the  angular 
momentum  transfer  j,  of  a  certain  Clebsch-Gordan 
coefficient  squared  multiplied  by  a  standard  angular 
function,  both  found  in  Table  I  of  I.  Assuming  the 
resonance  to  be  of  symmetry  A,  the  appropriate 
column  in  the  block  /  -  3  is  multiplied  by  the  angu¬ 
lar  functions  in  the  same  row  and  then  added  up  for 
al!  contributing  values  of  the  angular  momentum 
transfer,  j,  —  0,  2,  4,  and  6.  Then  the  following 
DCS's  are  obtained: 

for  fa{  A  =0)  , 

—r  =  0. 7 1 3<  1  4  5.69 cos:tf-  17.7 cos40 

ail 

(- 14. 3 cos6©)  ,  (la) 


for  fn(  A  =  I )  , 

~  0. 20 1 (  I  f  1 K.  2  cos’fl  -  52.1  cos40 
all 

+  40. 9 cosh0  )  ,  (lb) 
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ScoutMf^  Angle  I<J«g) 

FIG.  1.  Theoretical  angular  distributions  for  /-wave 
electrons  assuming  the  symmetry  for  the  molecular 
resonant  states. 

for  /&(  A  =2)  . 

=0. 534(  1+2.41  cos20  —13.6  cos40 
ail 

+  12.9cos60),  (lc) 

for  f<t>{  A  =  3)  , 

=  0.1I9<  I  +47.3 cos20 -  1 34 cos40 
all 

f99.9cos60).  (Id) 

In  Fig.  1,  these  four  curves  are  plotted  for  0=0° 
to  130°.  It  can  be  seen  that  they  are  distinctly  dif¬ 
ferent  from  each  other.  Generally  speaking,  the  two 
curves  fir  and  fi  are  too  high  at  60°  and  the  curve 
of  fh  it  too  low  at  40°  to  fit  the  experimental  data  on 
N.  and  CO  at  about  20  eV  (see  Figs.  2,  3,  and  5). 
However,  the  fa  curve  fits  some  data  very  well.  A 
remarkable  fit  shown  in  Fig.  2  is  the  N2  data  of 


i 


O 

P  V,  60  X  •?C 

Scotle'  I  nq  Angle  (deg) 

FIG.  2.  Experimental  angular  distribution  for  electrons 
having  excited  the  v  =3  state  in  N2,  compared  with  sim¬ 
plified  theory. 

Pavlovic  ei  al 4  for  e'  =  3  at  22.5  eV.  It  should  be 
pointed  out  that  the  fit  may  be  very  poor  at  other 
energies.  Fortunately,  such  variation  with  energy  is 
contradicted  by  the  more  recent  measurements  of 
Tanaka  et  al1  and  of  Tronc  et  al*  To  compare 
with  their  data  in  detail,  we  need  a  more  general  ex¬ 
pression  for  the  2Z„+  resonance  than  just  from  the  / 
wave. 

We  now  derive  the  general  DCS  for  a  reso¬ 
nance  in  a  homonuclear  molecule.  In  principle, 
higher  partial  waves,  e.g.,  /  =5,7, . . .  may  also  con¬ 
tribute  to  the  resonance,  but  here  they  can  be  exclud¬ 
ed  from  consideration  of  the  large  centrifugal  bar¬ 
rier,  as  /  =4,6, . . .  were  excluded  from  the  2IIg  reso¬ 
nance  at  a  lower  energy.  Then  the  ungerade  symme¬ 
try  requirement  leaves  us  with  just  the  p  wave  (/  =  1 ) 
as  well  as  the  /  wave  (/  =3).  This  admixture  of  two 
partial  waves  in  a  resonance  can  be  described  by  a 
single  parameter,  called  the  mixing  parameter  tan/?. 
In  fact,  precisely  this  situation  arose  in  II  of  this 
series.10  We  find 


—  (p'*-p  )=troy„  2  ( (I  +  cos2/?  )2(  10,y,0 1  1O)20(/,1 1 1  !0)  +  (  1  —  cos2/?  )2(  30,7,0  |  3O>20(./„33330) 


+  2  sin22/?(30,y,0  |  IO)20(y„  1 1330)]  . 


(2) 


In  Ht|.  (2),  a,,-,,  is  the  integrated  vibrational  excitation  cross  section  which  may  be  regarded  as  a  normaliza¬ 
tion  constant  here.  The  summation  over  the  angular  momentum  transfer  j,  is  restricted  to  even  values  by  the 
Clebsch-Ciordan  coefficients  (j tm  >  \j\tn  i )  and  the  parity  favoredness  ( 2  +  -  * 21  + )  of  the  transition.  The 
angular  function  0  is  again  given  in  Table  I  of  I.  Applying  the  usual  procedure  of  evaluating  the  DCS  to  Eq. 
(2 1.  we  obtain 
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" ) = av •„  [  0. 1 5(  1  +  cos2/3  )2(  1  +  2  cos2# ) 


+ 0. 1 78(  1  -  cos2/3  )2(  1  +  5 . 69  cos2#  -  1 7. 7  cos4#  +  1 4. 3  cos6# ) 


+O.183(2sin220)(l  +  l.O9cos2#)J  . 


(3) 


The  first  term  in  Eq.  (3)  [as  in  Eq.  (2)]  represents 
the  pure  pa  wave.  The  second  term  is  due  to  the  / a 
wave  and  therefore  is  identical  to  Eq.  (la).  Finally, 
the  last  term  is  the  usual  interference  term.  Over 
the  width  of  the  resonance,  P  is  expected  to  vary 
only  slowly  with  the  incident  energy  E.  In  principle, 
P  can  be  obtained  from  an  ab  initio  calculation. 
Here  we  choose  to  regard  it  as  a  parameter,  obtained 
by  fitting  Eq.  (3)  to  the  experimental  DCS. 

The  available  vibrational  excitation  data 
(v—0-*  1)  in  N:  consist  of  the  results  of  Tanaka 
et  at.1  at  18,  20,  22.5,  and  25  eV,  and  of  Tronc 
et  al*  at  18,  22.5,  and  25  eV.  All  the  DCS’s  show  a 
forward  peak,  a  minimum  at  90°,  and  some  flatten¬ 
ing  at  50°.  However,  they  differ  in  detail  by  as 
much  as  50%  even  at  the  same  energy.  The  error 
bars  given  by  Tronc8  in  Fig.  3  at  18  eV  indicate  their 
estimate  of  the  relative  error  is  ±8%,  while  their  ab¬ 
solute  error  is  stated  as  30%.  In  the  work  of  Tana¬ 
ka  et  al.  the  estimated  relative  maximum  error  was 
given  as  25%.  Thus  the  data  often  lie  outside  of 
each  other’s  error  bars  if  they  are  all  normalized  to 
unity  at  90°.  In  any  event,  the  scatter  in  the  experi¬ 
mental  DCS  appears  to  be  too  large  for  fitting  Eq. 
(3)  to  determine  the  mixing  parameter,  in  contrast  to 
the  1  8-eV  resonance  data  in  CO  discussed  in  II.  In¬ 
stead,  we  try  to  find  a  value  of  P  in  Eq.  (3)  which 
will  produce  a  theoretical  angular  distribution  which 
possesses  all  the  consistent  features  of  the  experi¬ 
mental  DCS  described  above. 

In  Fig.  3  we  show  a  curve  (solid  line)  obtained 
from  Eq.  (3)  with  casip—  —0.5,  i.e.,  /3  =  60°  or  the 
mixing  parameter  tan0=1.73.  Explicitly,  the  DCS 
is  then  given  by 


=0.71<70)(  1+4.02  cos2#  — 10. 1  cos4© 

aft 


+  8. 05  cos6#).  (4) 

Note  that  it  displays  a  forward  peak,  a  flattening  at 
50°,  and  a  dip  at  90°.  (Symmetry  about  90°  is  re¬ 
quired  by  theory  as  shown  in  II.)  For  comparison, 
only  some  experimental  DCS’s  are  shown,  selected 
on  the  basis  of  larger  integrated  cross  sections 
and/or  more  angular  positions.  They  are  normal¬ 
ized  to  theory,  Eq.  (4),  at  50°.  In  Fig.  3,  the  experi¬ 
mental  DCS's  are  now  seen  to  be  reasonably  con¬ 
sistent  with  each  other  and  with  theory.  Thus  the 


experimental  DCS’s  of  Tronc  et  al.  and  Tanaka 
et  al.  from  18  to  22.5  eV  for  vibrational  excitation 
in  N2  consistently  exhibit  the  shape  of  a  2Z„+  reso¬ 
nance  with  the  ratio  of  /-  to  p- wave  amplitude 
equaling 

(1+0.5) 

(1-0.5) 

Nevertheless,  the  sharp  disagreement  with  the  data 
of  Pavlovic  et  al.  remains  to  be  explained.  We  will 
show  how  high  resolution  discriminates  against 
some  rotational  components  and  leads  to  different 
measured  DCS's  from  those  with  lower  resolution. 

So  far,  we  have  neglected  the  rotational  structure 
of  the  molecules.  At  first  sight  that  seems  entirely 
appropriate  since  the  rotational  constant  of  N2  is 
5  =0.25  meV,  while  the  energy  resolutions  are  quot¬ 
ed  as  50  meV  in  the  works  of  Tronc  et  al.  and  Tana¬ 
ka  et  al.,  and  30  meV  in  Pavlovic  et  al.  However, 
the  predominantly  /  resonance  implies  sizable  rota¬ 
tional  excitation  cross  sections  of  the  type 
A j  —  ±4,  ±6  as  well  as  A j  —  ±2,  and  the  pure  elastic 


FIG.  3.  Theoretical  and  experimental  angular  distribu¬ 
tion  for  electrons  having  excited  the  t>  =  1  state  in  N2  (ar¬ 
bitrary  units). 
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A y  0.  In  general,  1  ho  additional  energy  loss  due  to 
a  rotational  braneh  A  j  is 

[27  A;  4-  ( Ay  )2  +  Ay  ]fl  . 

At  room  temperature,  y  =  10  is  a  typical  rotational 
state  with  significant  population;  for  this  state,  the 
energy  losses  for  Ay  =4  and  Ay  =6  are  25  and  40 
meV,  respectively.  Now  we  see  that  the  DCS  of 

do 

dii 


Puvlovic  el  at.  essentially  excludes  rotational  excita¬ 
tion  of  the  type  Ay  =  6,  and  partially  those  with 
Ay  =4,  while  those  of  Tronc  el  al.  and  Tanaka  el  al. 
partially  exclude  those  with  Ay  =  6. 

To  scrutinize  this  issue  of  rotational  excitation  in 
addition  to  vibrational  excitation,  we  need  the 
lotaiional-vibrational  DCS.  With  suitable  modifica¬ 
tions  for  the  22u+  resonance,  Eq.  (2)  of  I  becomes 


( /  V— yt>)  =  irovv  2  O' 0,7,0  | y'0)2[  ( 1  +  cos20  )2<  10,y,0 1  10)2©(y„  11110) 
i, 


+  ( 1  —  cos20)2(3O,/,O  |  30)2©(/,, 33330) 

+  2  sin220<  10,/,  0  j  10H  30,/, 0  30)©( / „  13130) 

+  2  sin22/?(  30,/, 0  j  10)2©(  /,,  1 1 330 ) 

—  2  sin220(  10,/,0  j  30)(  30,/,0 1  IO)0(/„13310)]  .  (5) 


Accepting  the  value  of  —0.5  for  cos2/3  found  earlier,  we  can  now  calculate  each  individual  rotational- 
vibrational  cross  section  in  N2  for  all  significantly  populated  levels  in  N2  at  the  molecular  beam  temperature. 
For  example,  the  Ehrhardt  group9  measured  their  temperature  to  be  500 1 30  K.  Rather  then  calculating  all  26 
or  so  DCS’s  for  each  populated  rotational  state  and  averaging  the  population,  we  calculate  the  DCS  for  each 
rotational  branch  for  the  median  rotational  level.  For  N2  at  500  K,  the  median  level  is  j  =  10  (population 
8.56%)  with  the  population  of  j  =0  to  9  totaling  42.9%,  and  the  population  of  y  =  1 1  to  oc  totaling  48.5%. 
We  use  the  high-/  approximation  developed  elsewhere,"  which  in  effect  neglects  terms  of  order  y~2  in  compar¬ 
ison  with  unity.  Therefore  it  is  accurate  to  about  1%.  The  DCS  for  each  rotational  branch  is  then  given  by 

Ay  =0)  =  j<v„[  ( 1  +cos2/3)2(0. 15+  l.O5cos20) 

ail 

+  ( 1  —  cos2(3  )2(0.09  +  2. 74  cos20  —  8. 82  cos40  +  7. 1 9  cos60 ) 


+  2  sin220(0. 87-0.58  cos20  +  1.61  cos40 )]  , 


do 

dil 


(Ay =  ±2)  = 


1 


[  ( 1  +cos20  )2(0.23  +0.08  cos20 ) 


+  (1  —  cos20)2(  0.  1O+O.77cos20  —  1 . 99  cos40  + 1 . 38  cos60 ) 


(6a) 


+  2  sin2j5(0. 1 4 + 0. 50  cos20  -  0. 3 1  cos 40 )  ]  , 


do 

dn 


(Ay  =  ±4)  = 


i 


[  ( 1  -  cos20  )2(0. 13—0.17  cos20  +0.06  cos40  +  0. 1 3  cos60 ) 


(6b) 


+  2  sin22/3(  0. 23  —  0. 20  cos20  —  0. 03  cos40 )  ]  , 


(6c) 
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~{^j=± 6)=-j-cv„  1±-  ( 1  —  cos2/?>2( 0.73 +Q. 55 cos20  —0.28 cos4#  40.02  cos60)/lO  .  (6d) 


The  results  of  Eq.  (6)  may  be  directly  verified  when 
the  Ehrhardt  group9  reruns  their  experiment  and 
meticulous  analysis  at  an  energy  of  20  eV. 

For  the  moment,  we  only  want  to  explore  the  ef¬ 
fect  on  the  measured  vibrational  DCS  when  some 
rotational  branches  are  inadvertently  excluded 
through  achieving  higher  resolution.  In  this  respect, 
the  most  important  branch  is  the  Ay  =  4  6  located  at 
an  energy  loss  of  40  meV  with  a  relative  cross  sec¬ 
tion 

o(Ay  =6)/aMlm 

of  10%.  The  next  important  branch  is  the  Ay  =  +4 
located  at  25  meV  with  a  relative  cross  section  of 
about  11%.  For  convenience  their  DCS,  given, 
respectively,  by  Eqs.  (6c)  and  (6d),  are  shown  in  Fig. 
4.  Since  they  comprise  only  large  angular  momen¬ 
tum  transfer  j,  in  Eq.  (5),  they  are  nearly  isotropic, 
as  classically,  electrons  undergoing  large  change  in 
angular  momentum  can  be  expected  to  “forget” 
their  incident  direction.  Relatively  unimportant  are 
the  corresponding  superelastic  branches  with  the 
same  DCS,  Ay  =  —  6  and  —4,  located  at  —  19  and 
-  15  meV  with  relative  cross  sections  of  4  and  5  %, 
respectively.  It  is  apparent  that  with  improving 
resolution  more  of  these  near  isotropic  DCS’s  will  be 
excluded  in  the  measurement.  Hence  the  experimen¬ 
tal  DCS  will  tend  to  be  too  high  in  the  forward 
direction  as  is  the  case  in  the  data  shown  in  Figs.  2 
and  3.  Further,  if  the  experimental  resolution  is 
varied  with  the  scattering  angle,  the  measured  DCS 
can  take  on  an  unexpected  shape.  Such  appeared  to 
be  the  case  in  the  work  of  Pavlovic  et  al.*  when  the 
angle  was  fixed,  and  the  incident  energy  was  varied 
in  the  measurement.  The  best  energy  resolution 


°  70  t.O  90  t.'0 


FIG.  4.  Theoretical  angular  distributions  for  /-wave 
electrons  with  the  selection  of  specific  rota’ional  branches. 


achieved  was  30  meV  at  some  angle,  but  at  other  an¬ 
gles  the  resolution  might  have  been  lower.  In  con¬ 
trast,  Tanaka  et  al.1  fixed  the  energy  while  varying 
the  angle,  as  did  Tronc  el  al*  In  both  latter  experi¬ 
ments,  the  resolution  was  50  meV,  which  appeared 
low  enough  to  include  most  of  the  Ay  =  4  6  rotation¬ 
al  branch.  Therefore  the  measured  DCS  of  Pavlovic 
et  al.  can  be  expected  to  be  quite  different  from  the 
latter  two  experimental  groups  and  from  theory. 

Next,  we  turn  our  attention  to  the  equivalent  reso¬ 
nance  in  CO.  DCS’s  have  been  measured  by  Tronc 
et  al.*  at  16,  19.5,  and  23  eV,  and  by  Chutjian  and 
Tanaka6  at  20  eV.  Figure  5  shows  that  these  DCS’s 
are  all  quite  consistent  with  each  other  and  fit  the 
theoretical  fa  expression  (dashed  line)  given  in  Eq. 
(la)  quite  well.  Once  more  we  can  conclude  that  e- 
CO  vibrational  excitation  is  dominated  by  a  2X  + 
resonance  at  about  20  eV  with  a  width  of  at  least  3.5 
eV.  In  principle,  a  2£+  resonance  may  contain 
s(l  =0),  y>(/  =  l),  d(/=2)as  well  as  /(  /  =  3)  com¬ 
ponents.  However,  it  can  be  argued  that  only  the  d 
component  needs  to  be  considered  as  only  it  is  dipo- 
larly  coupled  to  the  dominant  /  component. 

Then  the  DCS  for  the  vibrational  excitation  of 
CO  is  essentially  the  same  as  for  N2  given  in  Eq.  (2), 
except  for  the  replacement  of  the  p  wave  by  the  d 
wave.  For  the  sake  of  completeness,  it  is  given 
below: 


0  30  60  90  i 20 

ScoHertfiq  AnQie  (d*q) 


FIG.  5.  Theoretical  and  experimental  angular  distribu¬ 
tions  for  electrons  having  excited  the  r  I  state  in  CO. 
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~~  ( v'*~u  )  =  nau'e  2 1  <  1  +cos20)2<  20,7,0 1  2O)20(7„  22220) 

all  i 

Jt 

+  ( 1  —  cos2/3)2(30J,0  |  30)2©U,, 33330 ) 

+  2  sin22j8(  30,7,0 1  2O)20( j, ,22330 )]  .  (7) 


In  comparison  with  Eq.  (la),  a  better  fit  with  all  data  is  obtained  with  cos2/?=  —0.75  in  Eq.  (8),  i.e., 

( u'«— u ) = ott„(Q. 78 )( 1  +  3. 59  cos20  —  11.7  cos40  +9.98 cos60 )  ,  (8) 

all 

as  shown  in  Fig.  5.  Nevertheless,  significant  discrepancy  between  the  DOS’s  of  Tronc  et  al.  (at  19.5  eV)  and 
Chutjian  and  Tanaka  up  to  20%  at  20  eV  prevents  us  from  obtaining  a  value  for  cos2/3  more  accurate  than 
— 0.75±0.25.  Thus  the  2£  +  resonance  in  CO  is  also  dominated  by  the  /  wave  as  in  the  case  of  the  22,+  reso¬ 
nance  in  N2.  However,  the  admixture  is  lower  in  this  case  with  the  d- wave  amplitude  ratio  to  the  /  wave  given 
by 

(1-0.75)  __  0.25  _  1 
( 1  -f  0.75)  1.75  —  7 

In  anticipation  of  the  measurement  of  the  rotational-vibrational  DCS  in  CO  at  20  eV  by  the  Ehrhardt  group, 
we  present  the  individual  differentia)  cross  section  below.  The  general  expression  from  the  initial  state  ( j,v )  to 
the  final  state  O', o')  is  given  by 

~(j'v'—jv)=vav-v  2 0°./«0  I/O)2 
j, 


X  [  ( 1  +cos20)2(  20,7,0 1  2O)200, 22220) 

+  ( 1  -  cos2/3  )2<  30,7,0 1  3O)20( 33330) 

-f2sin22^(20,7,0|  20)(30,7,0|  30 )0(;„23230) 

+  2  sin22/3(  30,7,0  1  2O)20( 7,, 22330 ) 

-2sin22/3(20, 7, 0|  30)(30,7,0|  2O)0O„23320)]  .  (9) 


In  practice,  as  in  the  case  of  N2,  only  the  rotational  branches  can  be  separated  out.  With  the  same  high-7  ap¬ 
proximation,  the  DCS  for  each  branch  is  found  to  be 

—  ( A7  =  0 )  =  7  tv,  [  ( 1  +  cos2/5  )2(  0. 4 1  - 1 . 69  cos20  +  2. 49  cos40 ) 
ail 

+  ( 1  -  cos20)2(0. 10  +  2.73  cos20  -8.81  cos40 + 7. 20  cos60) 


+  2  sin22£(  9. 34  cos0  —  2. 75  cos30  +  4. 48  cos50 )  ]  , 


^<A7  =  ±1)=7<V, 


( 2  sin22/3 )( 0. 20 — 0. 54  cos0  —  0. 38  cos20  +  2. 62  cos}0 


(10a) 


+0.67  cos40  —  2. 57  cos50 )  , 


(10b) 
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^.A;  =  ±2  )  = 


[  ( 1  +  cos2/?)2(  0.22 -0.30  cos2©  4  0.35  cos *0) 


+  ( 1  -  cos2 P )2<  0. 1 0  +  0. 77  cos2©  -  ! .  99  cos4©  +  1.38  cos "© ) 


^  i  Ay'=  i3)=|ct„1! 


“ii 


+  2  sin:2/?(  -  0. 08  cos©  +1.11  cos '©  +  1  08  cos'© ) )  , 
( 2  sin:2 P )(  0. 1 8  -  0. 1 2  cos©  0. 30 cos 2©  *-  0. 1 9  cos'© 


+ 0. 32  cos4©  -  0. 28  cos '© )  , 

[  ( 1  -i-cos2/J)2(0. 1 1  +0.08 cos’©  +0.01  cos4©  I 
+  ( 1  - cos2/3 )’(  0. 14  —  0. 17  cos’©  +  0. 07  cos4©  +  0. 1 3  cos*© ) 


(10c  i 


(lOd) 


+  2sin22/?<  -0.02 cos©  +0.11  cos'© -0.04 cos'©)]  , 


(lOe) 


( 2  sin22 P )(  0. 9 1  -  1 . 1 6  cos©  +  0. 66  cos2©  -  0. 28  cos  '© 


+  0. 20  cos4©  -  0. 33  cos'© )  / 1 0  . 


non 


and 


da  ...  . ,  i 

— -(A;=±6)=7<7„.1> 


( 1  -  cos2/?  )2(  0. 73  +  0. 55  cos2©  +  0. 28  cos4©  +  0.02  cosh0 )  / 1 0  . 


For  the  purpose  of  comparison  with  the  anticipated 
data  of  the  Ehrhardt  group,  the  initial  rotational 
state  j  should  be  taken  to  be  the  median  level  at  the 
molecular  beam  temperature.  At  500  K  in  the  Jung 
et  a! ?  experiment,  two  levels  tie  as  the  median 
j  =  10  and  1 1  with  a  relative  population  of  6.35  and 
6.15%,  respectively.  (The  sum  from  j  =0  to  9  is 
42.8%,  and  the  sum  from  j  =  12  to  oc  is  44.7%.) 

In  conclusion,  we  have  applied  a  previous  theory 
to  elucidate  the  DCS  of  an  /  resonance  in  electron 
scattering  from  a  diatomic  molecule.  The  20-eV  res¬ 
onances  in  N2  and  CO  are  both  confirmed  to  be  of 
22  +  in  symmetry  with  a  dominant  /  component. 
Present  in  the  former  is  a  p  component  with  a  rela¬ 
tive  amplitude  of  0.33,  and  in  the  latter  a  d  com¬ 
ponent  with  a  relative  amplitude  of  0.14. 

Through  analysis  of  the  rotational  branches,  we 


(I0g) 

I - 

have  shown  that  spurious  DCS's  may  be  obtained 
through  too  high  an  energy  resolution.  Even  a  mod¬ 
est  resolution  of  30  meV  used  by  Pavlovic  el  al  * 
will  exclude  some  rotational  branches,  resulting  in  a 
distorted  vibrational  DCS. 

Finally,  we  have  presented  the  DCS's  of  all  the 
rotational  branches  in  the  vibrational  excitation  of 
N2  and  CO  in  the  2 2  +  resonance  region  of  about 
16—25  eV.  They  await  direct  confirmation  by  the 
new  experimental  technique  developed  by  the 
Ehrhardt  group. 

The  author  thanks  H.  Sakai  for  useful  discus¬ 
sions.  Helpful  correspondence  with  H.  Tanaka,  M. 
Tronc,  and  H.  Ehrhardt  are  gratefully  received. 
This  work  was  supported  bv  Air  Force  Contract  No. 
FI  9628-8  l-K-0007. 
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Abstract 


The  molecules  in  the  lower  part  of  atnosphere  are  in  thermal  equilibrium  with  the  local 
surrounding.  The  infrared  emission  resulted  by  their  thermal  excitation  can  ho  studied  to 
determine  their  temperature  and  concentration  profile  in  the  atmosphere.  We  used  a  cryo- 
genically  cooled  Fourier  spectrometer  mounted  on  a  bal loon-horne  platform  to  study  the 
infrared  emission  spectrum  of  atmosphere  for  this  purpose.  The  experiment,  which  took 
place  on  October  /,  1981,  at  To!  Ionian  AFB,  produced  t  he  analyxable  data  for  more  than  I  wo 

hours  at  altitude  of  2/000  28000  m  along  the  horizontal  and  the  vertical  line  of  sight. 

The  spectral  data  extended  from  550  cm'*  to  1000  cm'*  covering  the  CO2  15m  bands  (Ava  ■■  1) 
with  a  resolution  of  0.2  cm'*.  The  data  were  found  interesting  on  two  accounts:  (1)  a 
continuum  background  feature  was  superimposed  with  the  molecular  emission  feature;  (2)  the 
radiance  level  of  tne  CO2  15u  hands  observed  along  the  horizontal  line  of  sight  varied  by 
a  much  greater  degree  than  exnocted. 

Introdnct i on 

The  infrared  radiative  transfer  characteristics  in  the  atmosphere*  are  principally  con¬ 
trolled  by  two  parameters,  the  local  abundance  of  the  molecules  responsible  for  the 
emission  or  the  absorption  process,  and  the  local  temperature.  The  third  parameter,  the 
local  pressure,  plays  a  secondary  importance  to  the  problem.  In  the  lower  part  of  atmo¬ 
sphere  including  the  stratosphere,  the  molecules  are  generally  in  thermal  equilibrium  with 
the  local  surrounding.  The  radiative  transfer  characteristics  under  the  local  thermal 
equilibrium  condition  could  be  solved  in  a  straightforward  calculation,  if  those  parameters 
are  known.2  In  the  stratosphere,  some  of  these  parameters  are  noorly  known  at  present 
because  of  experimental  difficulties.  We  designed  an  experiment  to  collect  the  infrared 
radiative  transfer  data  in  the  stratosphere  by  observing  the  infrared  emission  along  the 
horizontal  line  of  sight.  A  cryogenical ly  cooled  Fourier  interferometer  spectrometer 
mounted  on  a  balloon-borne  platform  was  used  for  the  experiment.3  The  balloon  flight  took 
place  on  October  7,  1981,  at  Holloman  AFB,  New  Mexico,  and  obtained  the  spectral  data 
taken  at  the  balloon  ceiling  altitude  of  a,  28500  m  for  more  than  two  hours . 11  The  instru¬ 
mentation  part  of  our  experiment  is  detailed  in  another  paper  presented  at  this  meeting.5 
In  this  paper,  our  discussion  will  be  focused  on  the  spectral  data  obtained. 

In  the  past,  the  spectral  study  of  the  atmosphere  was  done  either  with  the  absorption 
scheme  using  the  sun  as  its  infrared  source,6  or  with  the  satellite-borne  spectrometer 
looking  down  toward  the  earth's  surface.7  The  absorption  scheme  ties  its  path  on  the  sun's 
position.  The  path  for  an  adequate  absorption  required  in  the  quantitative  study  of 
stratosphere  must  be  long,  because  few  molecules  are  available  there.  The  measurement  is 
then  limited  to  a  short  period  during  sunrise  or  sunset.  The  molecules  suspected  of  the 
diurnal  concentration  variation  cannot  be  studied  well  if  their  observation  requires  a  long 
absorption  path.  The  satellite-borne  experiment  measures  a  cumulative  effect  of  the 
absorption  or  the  emission  along  the  path.  It  requires  an  elaborate  analysis  scheme  if  the 
radiative  transfer  characteristics  of  a  thin  atmospheric  layer  is  sought.  The  scheme  that 
we  adopted  observes  the  atmospheric  emission  along  the  horizontal  line  of  sight  independent 
of  the  solar  position.  The  radiative  transfer  characteristic  of  a  thin  atmospheric  layer 
is  more  readily  available  in  our  measurement  scheme.  The  observation  extended  a  long 
period  of  time  over  the  entire  balloon  flight,  since  the  line  of  sight  for  observation  did 
not  tie  into  the  solar  position.  The  radiance  level  was  very  low.  At  the  maximum,  it 
would  correspond  to  the  blackiody  radiance  at  the  stratosphere  local  temperature  of  220°  ■v 
230°K.  The  sensitivity  of  measurement  was  a  critical  issue  in  our  experiment  and  it  is 
fully  discussed  in  our  other  paper. 


Balloon  flight 


The  flight  profile,  given  in  terms  of  the  balloon  altitude  as  a  function  of  time,  was 
monitored  using  radar  on  the  ground.  It  is  shown  in  Fig.  1  with  the  local  temperature 
measured  with  a  thermistor  bo  ometer  hung  below  the  gondola.  The  measured  value  of  the 
total  radiance  level  is  also  >h>wn  in  the  figure.  Even  though  the  laser  signal  failure 
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which  occurred  during  some  part  of  the  flight  resulted  in  some  loss  of  the  data,  the  inter¬ 
ferometer  produced  the  analyzable  data  for  more  than  two  hours  at  the  maximum  balloon 
altitude.  The  spectral  data  obtained  during  the  flight  consisted  of  three  parts:  (1)  the 
data  along  the  horizontal  line  of  sight;  (2)  the  data  along  the  vertical  line  of  sight;  and 
(3)  the  radiance  calibration  spectra  obtained  using  an  on-board  blackbody  source  which  was 
presumably  in  thermal  equilibrium  with  the  ambient  As  seen  in  the  thermistor  bolometer 
data,  the  ambient  temperature  at  the  altitude  of  28000  m  -v,  28500  m  fluctuated  between  270°K 
and  250°K,  closely  correlating  with  the  balloon  motion.  The  data  indicated  that  the 
stratospheric  ambient  temperature  measured  with  the  on-board  thermistor  was  considerably 
higher  than  that  of  the  model  atmosphere,  and  that  its  fluctuation  was  unexpectedly  large. 

Our  balloon  flight  was  accompanied  by  two  radio-sonde  flights.  The  temperature  measured 
by  the  radio-sonde  sensor  was  at  most  230°K,as  seen  in  Fig.  2,  at  this  altitude  range,  much 
lower  than  the  temperature  measured  with  the  on-board  thermistor.  Even  if  a  rapid  motion 
of  the  radio-sonde  balloon  was  accounted  for,  the  discrepancy  among  those  three  temperature 
readings  defied  a  single  explanation.  This  puzzle  remains  unsolved  at  present.  Both  data 
showed  the  temperature  fluctuation  of  approximately  10°K  at  the  altitude  of  28000  m.  The 
variation  in  the  total  radiance  level  measured  along  the  horizontal  line  of  sight  was  also 
found  an  unexpectedly  large  value  of  50%.  It  was  also  found  to  be  correlated  with  the 
measured  temperature. 

Calibration  of  the  radiance  level 

Two  sets  of  spectra  were  taken  using  the  on-board  blackbody  source,  one  at  17:30  GMT  and 
another  between  18:45  and  19:30  GMT.  Both  sets  were  taken  at  the  balloon  ceiling  altitude 
of  -v  28000  m.  Their  spectral  characteristics  as  well  as  their  radiance  levels  were  found 
very  similar.  We  estimated  these  spectra  to  be  those  of  a  255°K  blackbody  radiance,  and 
proceeded  to  calibrate  all  the  measured  spectral  data,  horizontal -looking  and  down-looking, 
accordingly.  The  calibrated  down- looking  spectra  showed  the  radiance  level  at  the  center 
of  the  CO2  15p  band  corresponding  to  215°K  ■c  220°K  and  that  of  the  window  region  to  the 
ground  temperature  of  285°K.  The  radiance  level  obtained  for  the  down-looking  spectra  was, 
therefore,  consistent  with  the  atmospheric  radiative  transfer  expected  for  this  spectral 
region . 

The  radiance  level  of  the  interferogram  signal  was  estimated  with  respect  to  the  noise 
level  contained  in  it.  Since  the  gain  factor  of  the  detector  amplifier  was  not  transmitted 
to  the  ground,  it  had  to  be  determined  from  the  recorded  interferogram  data.  The  noise 
level  in  the  interferogram  data  increased  or  decreased  whenever  the  detector  amplifier 
switched  its  gain  setting.  Since  we  had  no  reason  to  believe  that  the  absolute  level  of 
noise  in  the  interferogram  signal  varied  from  time  to  time,  we  determined  the  gain  of  the 
detector  amplifier  using  the  noise  level  observed  in  the  recorded  interferogram  signal .  We 
observed  the  noise  level  stepped  in  accordance  with  the  step  in  the  amplifier  gain  setting. 
Therefore,  our  procedure  for  estimating  the  detector  amplifier  gain  was  self-consistent. 

The  total  radiance  level  received  by  the  spectrometer  was  proportional  to  the  central  maxi¬ 
mum  modulation  and  it  was  estimated  by  multiplying  the  central  maximum  modulation  with  the 

fain  factor  which  was  determined  from  the  noise  level  observed  in  the  interferogram  very 
ar  from  the  central  maximum. 

An  exception  to  the  procedure  for  the  radiance  calibration  had  to  be  observed  for  the 
data  taken  after  19:30  GMT,  when  the  blackbody  source  was  removed  from  the  field  of  view  at 
the  last  time  in  the  flight.  The  radiance  level  determined  using  the  procedure  above 
showed  a  dramatic  increase  for  the  horizontal-looking  spectra  taken  after  19:30  GMT,  while 
the  balloon  altitude  and  the  measured  ambient  temperature  remained  unchanged  as  before.  At 
that  time,  the  liquid  helium  level  in  the  detector  dewar  became  definitely  low  and  the 
detector  temperature  started  to  rise.  We  reasoned  that  the  apparent  increase  of  the 
radiance  level  after  19:30  GMT,  shown  by  a  solid  curve,  was  resulted  by  an  improved  noise 
figure  of  the  detector,  rather  than  by  an  increase  in  the  incident  radiance  level.  The 
radiance  level  probably  remained  at  the  same  level  as  indicated  by  a  dotted  curve.  Anyway, 
the  data  obtained  after  19:30  GMT  certainly  had  a  larger  degree  of  ambiguity  than  the  rest 
because  the  radiance  level  calibration  was  not  repeated  after  the  observation  of  these 
spectra.  The  rest  of  our  discussion  will  be  made  excluding  the  data  collected  after  19:30 
GMT. 

The  sensitivity  of  our  spectrometer  was  determined  using  the  blackbody  calibration 
spectrum.  As  described  above,  the  calibration  source  was  assumed  at  255°K.  The  obtained 
blackbody  spectrum  and  the  determined  radiance  level  given  as  a  function  of  wavenumber  are 
shown  in  Fig.  3.  The  obtained  curves  were  used  to  determine  a  relative  radiance  level  of 
the  data  at  each  wavenumber  position.  The  absolute  value  of  the  radiance  level  was  then 
determined  by  normalizing  the  integrated  relative  radiance  value  to  the  total  radiance 
level  obtained  using  the  procedure  described  above. 


C02  bands 


The  major  atmospheric  bands  observable  in  our  measurement  are  CO2  (Avo  =  1)  transitions. 
All  bands  have  a  strong  Q  branch.  Table  I  lists  the  band  center  for  the  observable 
transitions.  Of  these,  the  Q  branch  of  the  (01101-00001)  band  of  the  626  isotope  should 
show  a  strong  saturation  even  at  our  balloon  altitude  with  an  assumption  that  the  COy 
concentration  in  the  stratosphere  is  given  by  the  mixing  ratio  of  3.3  x  10"5.  Our  obser¬ 
vation  along  the  horizontal  line  of  sight  was  made  at  the  elevation  angle  of  7.5°.  With 
the  CO2  mixing  ratio  of  3.3  x  10"^  and  the  mid-latitude  model  atmosphere,  we  can  expect 
the  CO2  column  density  to  be  approximately  5.0  x  10^0  molecules/cm2 . 


Table  I. 

Major  COo  Bands  in  the 

700  cm  *  Region 

Band 

Isotope 

Band  Center 

01101-00001 

626 

667.4  cm"* 

10001-01101 

626 

720.8 

10002-01101 

626 

619.8 

02201-01101 

626 

667.8 

03301-02201 

626 

668.1 

01101-00001 

636 

648.5 

01101-00001 

628 

662.4 

The  down-lookinR  data 

The  emission  spectra  observed  along  the  vertical  line  of  sight  is  presented  in  Fig.  A. 
All  data  observed  along  this  direction  varied  very  little  among  them  in  their  spectral 
characteristics  and  their  radiance  level.  As  mentioned  earlier,  the  source  radiance  level 
in  the  675  cm" *  region  was  found  to  be  t  4.0  x  10"°  W/cm^  str  cm'*,  the  blackbody  radiance 
level  of  210°  ^  220°K.  The  level  in  the  800  cm'*  region  was  1.1  x  10"*  W/cm2  str  cm"*, 
consistent  with  the  concurrently  measured  ground  temperature  of  285°K.  The  theoretical 
spectrum  calculated  using  the  FASC0D1B  code2  with  the  AFGL  atmospheric  line  parameters8  is 
shown  in  Fig.  5.  No  significant  differences  were  found  between  the  observed  and  the 
calculated. 

The  horizontal- looking  data 

The  data  observed  along  the  horizontal  line  of  sight  showed  their  spectral  characteris¬ 
tics  almost  identical  to  each  other.  Their  difference  was  noticed  in  the  radiance  level 
as  demonstrated  in  the  data  of  Figs.  6  and  7,  which  were  taken  at  18:15  and  18:45  GMT. 

Two  distinctive  features  were  found  in  these  observed  spectra:  (1)  a  continuum  feature 
which  amounted  to  the  radiance  level  of  2  '  3  x  10"°  watt/cm2  str  cm"*  at  the  600  cm"* 
region,  and  (2)  a  large  fluctuation  in  the  observed  radiance  level. 

A  theoretical  spectrum  shown  in  Fig.  8  was  generated  for  the  emissivity  of  those  C02 
transitions  observable  in  the  spectral  region.  The  calculation  assumed  that  a  column  of 
U  x  10zo  CO2  molecules/cm2  was  at  a  kinetic  temperature  of  230°K.  If  the  temperature  was 
different  from  230°K  by  a  small  degree,  the  emissivity  would  have  changed  little.  Thus 
the  curve  is  a  very  convenient  guide  for  estimating  the  radiance  level.  In  the  figure,  a 
solid  line  representing  a  radiance  level  of  5  x  10"°  watt/cmz  str  cm"*  gives  the  scaling 
factor  for  230°K.  If  the  temperature  is  210°K,  the  radiance  level  would  have  to  be  scaled 
by  a  dotted  line  which  provides  a  radiance  level  of  4.0  x  10"°  watt/cm2  str  cm'*.  Our 
observation  at  altitude  of  28000  m  with  the  elevation  angle  of  7.5°  should  have  yielded 
the  CO2  emission  with  its  column  density  of  5  x  102**  molecules/cm2.  The  Q  branch  of  the 
(01101-00001,  626)  band  was  expected  to  show  a  complete  saturation  (emissivity  =  l)  at  its 
peak.  The  radiance  level  at  this  peak  was  found  to  be  3.7  x  10"°  watt/cm2  str  cm"*  for 
the  spectrum  of  Fig.  7,  and  3.0  x  10"°  watt/cm2  str  cm"*  for  Fig.  6.  The  spectrum  of 
Fig.  7  could  be  interpreted  as  the  emission  at  210°K.  This  interpretation  ran  into  a 
difficulty  because  it  produced  no  satisfactory  explanation  for  the  continuum  feature.  The 
radiance  level  of  3.0  x  10"°  watt/cm2  str  cm"*  compounded  our  difficulty  in  explaining 
these  results. 

The  theoretical  calculation  obtained  using  the  FASC0D1B  code  with  the  CO2  column  density 
of  1.2  x  102*  molecules/cm2  at  220°K  is  shown  in  Fig.  9.  Both  theoretical  results  con¬ 
tained  no  indication  of  the  continuum  feature.  We  engaged  in  an  extensive  search  for  any 
experimental  artifacts  which  could  possibly  produce  a  continuum  feature,  and  found  none. 

We  are  rather  confident  that  the  continuum  feature  seen  in  Figs.  6  and  7  is  real  and  of 
atmospheric  origin,  even  though  we  have  no  explanation  for  it  at  this  moment. 

The  two  unsolved  mysteries  are  our  principal  discoveries  of  our  latest  balloon  flight. 


We  are  honing  that  future  flights  will  bring  some  light  onto  the  observed  mysteries  and  that 
we  will  gain  more  insight  into  the  stratosphere. 
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Figure  1.  Flight  Profile:  Balloon  Altitude,  Ambient  Temperature, 

and  Observed  Total  Radiance. 
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Abstract.  Analysis  of  rotational  branch  structure  in  high  resolution  electron  energy-loss 
spectroscopy  is  greatly  facilitated  by  a  high  /  approximation  to  the  Clebsch-Gordan 
coefficient.  A  more  accurate  approximation  than  Read's  is  found  for  molecules  in  the  I 
state.  When  assimilated  into  the  theory  of  Chang,  the  resulting  expressions  for  the  branch 
differential  cross  sections  for  N2  and  CO  turn  out  to  be  easily  evaluated  by  hand.  In  both 
cases,  good  agreement  is  found  with  experiment  in  both  the  shape  and  the  magnitude  of 
the  cross  sections  at  an  electron  energy  of  around  2  eV. 


With  the  continuing  improvement  of  energy  resolution  in  electron  energy-loss  spectros¬ 
copy,  it  has  finally  become  possible  to  separate  out  the  rotational  branches  in  typical 
molecules  such  as  N2  and  CO.  The  process  of  deconvoluting  the  branches  was  greatly 
facilitated  by  utilising  the  high  /  (rotational  quantum  number)  approximation,  elegantly 
advanced  by  Read  (1972).  However,  the  approximation  is  good  only  to  order  /"’,  or 
alternatively  by  Read’s  own  estimation,  the  error  is  within  a  factor  of  two  given  by 
(50/m„//)%.  Either  way  yields  an  error  of  10%  in  the  present  case  where  j  - 10. 
Further,  even  with  the  above  approximation.  Read's  expression  for  the  differential 
cross  section  is  sufficiently  cumbersome  as  to  require  a  computer  for  its  evaluation. 

In  this  letter,  an  alternative  high  j  approximation  is  proposed,  which  is  valid  to 
order  j~2  and  therefore  1%.  When  incorporated  into  the  theory  of  Chang  (1977a, 
b,  referred  to  as  I  and  II  respectively),  it  produces  simple  expressions  for  the  rotational 
branches  of  the  differential  cross  section,  which  are  easily  evaluated  by  hand.  To 
demonstrate  the  utility  of  the  present  method,  we  evaluate  these  expressions  for  N2 
and  CO  in  the  n  resonance  region  around  2  eV.  Comparison  will  be  made  with  the 
experimental  data  of  Jung  et  al  (1982)  and  the  theory  of  Read  and  Andrick  (1971). 
To  emphasise  that  the  present  theory  predicts  more  than  just  the  angular  distribution, 
comparison  will  also  be  made  for  the  integrated  cross  sections. 

In  contrast  to  the  theory  of  Read  (Read  and  Andrick  1971,  Read  1972),  the  theory 
of  Chang  contains  only  one  squared  Clebsch-Gordan  coefficient  which  depends  on 
the  rotation  quantum  numbers  /  and  /'=/'  +  A.  For  2  states,  all  magnetic  quantum 
numbers  vanish.  Thus,  we  only  need  a  high  /  approximation  for  the  following 

(/0,fO'/  +  AO)2  =  (2/  +  2A+l)('  '  '*A)2.  (1) 
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A  well  known  formula  (Edmonds  19S7)  gives  the  expression  in  terms  of  factorials. 


(/0,  rO|/  +  AO)2  = 


(2/  +  2A  + 1)(2;  +  A  -  /)!  (f  +  A)!  (f  -  A)! 

1  i 

(/  + 

' 

(2/+A+f  +  l)» 

+ 

> 
N>  | 

1*1 

"1 

(2) 


After  some  manipulation  we  obtain 
I/O,  tO\j  +  AO)2  =  ( 1 + 


u+A.1)!!U.A„1)!!(l+^)(1+-zl)...(1+--) 


2 i 


2t-)!(-)!(l+^)(1+-lzl)...(1+-tl)- 

(3) 

In  the  high  j  limit,  we  may  keep  terms  only  to  order  /”'  and  ignore  terms  of  0(/~2). 
Then  we  find 

/  A\  It—  A_1\fW,J-A_  mi 

(4) 


A-  l)!!(r  + A- 

1)H 

1  * 

Ml 

fr) 

1  • 

Note  that  the  dependence  on  j  in  equation  (4)  entirely  lies  in  the  trivial  factor  (1  +  A/2/). 
Further,  detailed  balance  is  satisfied  to  the  same  order  since 

To  confirm  the  validity  of  our  approximation,  the  right  hand  side  of  equation  (4) 
is  evaluated  numerically  for  /  =  10  in  table  1.  The  exact  values  below  the  approximate 

Table  1.  High  /  and  exact  squared  Clebsch-Gordan  coefficients  for  j  =  10. 
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0(-n 

1 
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/0.179N 

0.141 
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\0.18l) 

0.144 

1  (-2) 

0.525 

/0.338N 
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/0.141\ 
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\0.338/ 

0.198 

'0.144/ 

2  (-3) 

0.413 

/0.265\ 
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\0.265/ 

0.173 

3  <-41 

0.359 

/0.219\ 
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ones  are  obtained  from  the  tables  of  Rotenberg  et  al  (1959).  The  agreement  is  about 
2%,  which  is  hardly  surprising  since  we  have  neglected  terms  of  order  /2  with  /  =  10. 
The  sum  rule,  which  requires  that  the  columns  add  to  unity,  is  satisfied  by  the 
approximate  values  (as  by  the  exact  ones)  to  ±0.001.  Finally  we  note  that  the  double 
factorial  in  equation  (4)  is  related  to  the  gamma  function  by  (Abramowitz  and  Stegun 
1964) 

(2n  - 1)!!  =  1.3  . . .  (In  - 1)  =  2nr (n  +  2)/v^.  (6) 


So  (-1)!!  has  the  value  of  unity  when  n  =  0  in  equation  (6). 

Now  the  approximation  in  equation  (4)  will  be  utilised  in  the  analysis  of  branch 
structure  in  N2  and  CO  at  the  2  eV  resonance.  For  the  n  resonance,  the  expression 
for  the  rotational-vibrational  differential  cross  section  has  been  given  in  equation 
(II.4).  In  the  case  of  N2,  the  2ng  resonance  implies  that  cos  20  =  -1,  resulting  in  this 
simple  expression, 


dcr  . 

s("'- 


■  v'j  +  A)  =  4tt<tl 


-K)? 


(t-A-  l)i!(f  +  A- 1)!! 


(21,  rO|21)20(r;  22220). 


(7) 


In  equation  (7)  crm ■  may  be  regarded  as  a  normalisation  factor,  which  is  constant 
except  for  the  often  negligible  factor  of  [(E  -£a«>.4/)/£]2  5  for  d  waves.  It  is  obvious 
that  the  dependence  on  /  lies  only  in  the  first  bracket  and  so  members  of  the  same 
rotational  branch  A  have  the  same  angular  distribution.  For  the  analysis  of  energy- 
loss  data,  Jung  et  al  (1982)  had  to  represent  the  entire  branch  by  a  single  ‘level' 
labelled  by  the  appropriate  typical  value  of  /.  Evidently,  this  value  should  be  the 
median  so  that  about  half  of  the  rotational  population  lay  on  one  side  and  half  on 
the  other.  In  terms  of  the  area  under  the  energy-loss  spectra,  the  same  choice  should 
be  made  with  a  slight  shading  to  the  low  j  side  for  A  >  0  (and  the  high  /  side  for  A  <  0) 
because  of  the  (1  +  A/2/)  factor.  At  a  temperature  of  500  K  in  the  Jung  experiment, 
the  median  value  of  j  for  N2  is  10  with  42.9%  of  the  population  in  /  =  0  to  9  and 
48.5%  in  J  =  1 1  to  oo. 

We  proceed  to  evaluate  equation  (7)  explicitly  for  the  case  of  N2.  The  Clebsch- 
Gordan  coefficients  and  the  angular  functions  0  may  be  found  in  table  1  of  I.  The 
results  are 


T?(ty  -» »7)  =  <Tw  (0.301 

all 

-  1. 540  cos2  0  +  2.3 10  cos4  0) 

(8  a) 

~(ty  -» v'j  ±2)  = 

1  ±^(0.651  -0.663  cos2  0  +0.886  cos4  0) 

(8b) 

^(ty  -»  v'j  ±  4)  = 

1  ±-)(0.503  +0.336  cos2  0  +  0.055  cos4  0). 

(8c) 

After  summing  over  r.ll  rotational  branches,  we  obtain 


— (v  -*t>')  =ow(0.536- 1.607  cos2  0  +  2.5  cos4  0)  =  f|crTO(  1  -  3  cos2  0  +  ^  cos4  6) 


(8d) 


which  first  appeared  in  the  work  of  Ehrhardt  et  al  ( 1 968).  For  the  vibrational  excitation 
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v  =  0-»  v'  =  1,  say  at  2.27  eV,  our  results  based  on  equations  (8)  are  indistinguishable 
from  those  of  Read  and  Andrick  (1971),  where  Read's  high  j  approximation  has  not 
been  used.  As  discussed  by  Jung  et  al  (1982),  experiment  is  in  excellent  agreement 
with  theory.  However,  one  minor  difference  is  the  ratio  of  areas  under  the  rotational 
branches.  Our  analysis  gives  the  ratios  for  A  =  -  2  to  +2  as  0.82  and  for  A  =  -4  to 
+4  as  0.63,  while  that  of  Jung  et  al  obtained  from  Read's  approximation  gives  the 
values  0.76  and  0.56.  Further,  it  should  be  emphasised  that  equations  (8)  contain 
only  one  normalisation  factor  a0i,  so  that  once  it  is  determined  from  any  one  of  the 
four  expressions,  the  other  three  show  excellent  agreement  with  experiment,  not  only 
in  the  shape  but  also  in  the  magnitude  of  the  differential  cross  sections. 

The  case  of  the  2fl  resonance  has  already  been  treated  by  Chang  (1977b)  in  II. 
With  the  high  /  approximation  developed  here,  we  obtain  the  differential  cross 
sections, 

~^(t >j  -  v'j)  =  (ow)[(l  +cos  2/3)2(0.038  +  1.103  cos2  6) 
ail 


+  ( 1  -  cos  2/3)2(0.298  - 1.532  cos2  0  +  2.302  cos4  0) 
+  (2  sin2  2/3)(-0.507  cos  0  + 1.478  cos3  0)] 


do-  . 

dn(”' 


■v'i±\) 


(9a) 


=  (ow)^l  ±^2  sin2  20 

x  (0. 102  +  0.319  cos  0  +  0.230  cos2  0  -  0.650  cos3  0) 

%(vj^v'i±2) 

=  (aw/ 10) (l  [U  +cos  20)2(O.56+O.19  cos2  0) 

+  ( 1  -  cos  20  )2(0.67  -  0.66  cos2  0+0.89  cos4  0 ) 

+  2  sin2  20(-O.14  cos  0-0.40  cos3  0)] 
da 

—  (e/-*e/±  3) 


(9b) 


(9c) 


=  (al„,/10)(l±|r)2sin2  20 

x  (0.54  -  0.67  cos  0  +  0.27  cos2  0  -0. 1 3  cos3  0)  (9 d) 

^(u/-*u'/±4)  =  (a„u/10)^l±^(l-cos  20)J(O.5O  +  O.33  cos2  0+0.06  cos4  0). 

(9c) 

After  summing  over  all  the  branches,  we  obtained  equation  (II. 2)  as  expected. 

Figure  1  shows  the  fit  of  equation  (II.2)  to  the  vibrational  excitation  data  of  Jung 
et  al  (1982).  They  obtained  a  value  of  cos  20  =0.15,  which  differed  somewhat  from 
the  previous  value  of  0.25  resulting  from  fitting  to  the  data  of  Ehrhardt  et  al  (1968). 


Letter  to  the  Editor 


L877 


Scattering  onglf  «3»g! 


Flfwcl.  Differential  crocs  sections  by  rotational 
branches  in  CO  of  resonant  vibrational  excitation 
at  1.8  eV.  Experimental  points  are  from  Jung elal 
(1982),  while  theory  is  from  equations  (9)  and 
(11.2),  with  cos  20  =  0.15. 


Also  seen  are  the  differential  cross  sections  for  the  rotational  branches  A  =  0,  and 
+1.  The  agreement  between  theory  and  experiment  is  better  revealed  here  than  in 
the  ratios  shown  by  Jung  et  al.  The  other  branches  given  by  equations  (9c  to  e)  are 
virtually  indistinguishable  from  the  data  and  are  not  shown.  Nevertheless,  significant 
discrepancy  exists  in  the  forward  direction  in  all  branches  of  CO,  which  is  even 
apparent  in  the  case  of  N2.  Examination  of  the  energy-loss  profiles  of  Jung  et  al 
(1982)  at  small  angles  shows  that  virtually  all  broadening  is  due  to  the  instrumental 
width  rather  than  the  inelastic  branches.  Therefore,  the  deconvolution  technique 
become  unreliable  and,  in  addition,  some  discrepancy  may  also  be  due  to  the  non¬ 
resonant  contributions  neglected  by  theory  as  discussed  by  Jung  et  al  (1982).  In 
contrast  to  our  theory,  Jung  et  al  found  that  the  theory  of  Read  and  Andrick  (1971) 
did  not  agree  well  with  their  data,  in  spite  of  the  presence  of  more  parameters. 

Finally,  our  theory  is  capable  of  predicting  the  integrated  branch  ratios  to  the 
vibrational  cross  section,  i.e.,  <r(i>/  -*v'j  +  A)/o-(t)  -*v').  Table  2  shows  these  ratios 
from  theory  and  from  the  experiment  of  Jung  etal  (1982).  In  N2,  because  all  differential 
cross  sections  are  symmetrical  about  90°,  they  are  able  to  account  reliably  for  the 
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Table  2.  Integrated  rotational  branch  ratiost  in  N2  and  CO  for  vibrational  excitation 


Theoryt  Experiments 


Theory}  Experiment!! 


t  Ratios  are  taken  with  respect  to  the  vibrational  excitation  cross  section  cr(w  =  0  -*  1 ). 
Each  column  should  add  to  1 . 

t  Strictly  speaking,  each  theoretical  ratio  should  be  multiplied  by  the  flux  conservation 
factor  [(E-  E„i)/(E  -E„0)1'2"ii/2,  which  differs  from  unity  by  less  than  3%  in  the  present 
cases. 

S  Jung  el  al  (1982)  2.2S  eV:  for  conversion  into  absolute  cross  sections,  the  denominator 
is  5.4 1  x  10  IA  cm2.  Ratios  at  other  resonant  energies  are  virtually  identical. 

||  Jung  el  al  (1982)  1.8eV,  unavailable,  (missing  data  from  120  to  180°):  for  conversion 
into  absolute  units,  the  denominator  is  7.4S  x  10" 16  cm2. 


missing  range  of  120  to  180°.  It  can  be  seen  that  the  agreement  between  theory  and 
experiment  is  about  10%.  For  the  case  of  CO,  the  differential  cross  section  of  the 
rotational  branches  are  quite  asymmetrical  about  90°,  so  no  experimental  values  are 
given.  Nevertheless,  the  vibrational  (summed  over  rotational  branches)  cross  section 
is  symmetrical,  so  a- (0-»  1)  at  1.8  eV  can  be  determined  reliably  to  be  7.45  x  10_I6cm2. 
(This  value  is  about  twice  as  large  as  the  one  measured  by  Ehrhardt  et  al  (1968)).  So 
theoretical  absolute  rotational  cross  sections  for  each  branch  can  be  obtained  by 
multiplying  this  value  to  the  ratios  given  in  the  last  column  of  table  2. 

In  summary,  the  high  /  approximation  developed  here  is  shown  to  be  more  accurate 
than  Read’s  by  one  order  of  magnitude.  When  incorporated  into  the  theory  of  Chang, 
the  resulting  expressions  are  so  simple  as  to  obviate  the  use  of  a  computer. 
Theoretical  rotational  branch  cross  sections  in  vibrational  excitation  in  N2  and  CO 
are  evaluated  and  shown  to  be  in  good  agreement  with  experiment  in  both  shape  and 
magnitude.  Application  of  theory  to  non-£  states  is  deferred  to  a  future  work. 

The  author  acknowledges  useful  correspondence  with  H  Ehrhardt,  and  thanks  J  F 
Walker  for  a  fruitful  discussion.  This  work  was  supported  by  Air  Force  Contract 
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Abstract.  We  have  investigated  the  vacuum  ultraviolet  absorption  of  N2  at  77  K.  The 
bands  at  810,  81$  and  825  A  are  identified  as  the  (0,0)  bands  of  the  7f,  6f  and  Sf 
complexes,  each  with  O,  Q  and  S  branches.  Theoretical  analysis  of  our  data  confirms  a 
theoretical  value  for  the  quadrupole  moment  of  2.1  e  a and  predicts  a  value  of  19  a,', 
for  the  polarisability  of  the  N  \  (v  =  0)  core. 


Typical  molecular  Rydberg  states,  which  are  core  penetrating  are  usually  studied  in 
absorptions  from  the  ground  state  in  the  wavelength  region  below  1000  A  (Yoshino 
el  al  1976,  1979).  On  the  other  hand,  non-penetrating  Rydberg  states  are  found  in 
infrared  emissions:  5g-4f  in  NO  by  Dressier  el  al  (1981)  and  5g-4f  and  4f-3d  in  H2 
and  D2  by  Herzberg  and  Jungen  (1982).  We  report  here  absorptions  from  the  ground 
state  to  the  non-penetrating  nf  Rydberg  complexes  of  N2.  The  spectra  for  each 
complex  resemble  the  P,  Q  and  R  branches  of  an  ordinary  electronic  band,  except 
that  the  spacing  between  adjacent  lines  here  is  4 B  rather  than  2 B  (where  the  rotational 
constant  B  is  nearly  equal  for  both  the  upper  and  lower  states).  Our  theory  attributes 
this  feature  to  all  homonuclear  molecules  in  a  £  ground  state.  Therefore  it  was  not 
present  in  the  4f  and  5f  absorption  study  of  the  2n  state  of  NO  (Jungen  and  Miescher 
1969),  which  was  complicated  further  by  the  unpaired  electron  spin. 

In  this  work  the  absorption  spectra  of  nitrogen  were  photographed  by  a  6.65  m 
McPherson  Model  265  vacuum  spectrograph  in  the  second  order.  The  reciprocal 
linear  dispersion  of  the  spectrum  is  approximately  0.63  A  mm"’.  Figure  1  (plate) 
shows  the  absorption  spectrum  of  N„<  at  liquid  nitrogen  temperature  from  A  = 
808-830  A.  Features  marked  nf  clearly  correspond  to  the  approximate  wavelengths 
of  n=7,  6  and  5  in  (/-Ryd/n2),  if  we  adopt  the  ionisation  potential  of  /  = 
125  667.5  cm'1  (Loftus  and  Krupenie  1977).  Not  shown  is  the  appearance  of  the 
isotope  bands,  l5N2,  at  about  the  same  wavelength  which  confirms  the  vibration 
assignment  of  these  bands  as  the  (0, 0).  For  greater  clarity,  the  810  A  feature  is  shown 
in  figure  2(a)  on  an  expanded  wavenumber  scale.  On  the  left  side  of  this  figure,  it  is 
evident  that  the  main  absorption  lines  are  regularly  spaced  at  about  8  cm’1  apart, 
which  is  4 B  since  the  rotational  constants  of  N2  and  N2  are  B  =  1.9982  and  B'  = 
1 .93 1 9  cm  1  respectively.  In  fact,  the  intensities  are  seen  to  alternate  with  strong/weak 
lines  as  expected  in  an  ordinary  P  (or  R)  branch.  However  as  just  noted,  the  spacing 
here  is  4 B  rather  than  2 B,  so  by  analogy  we  label  this  group  the  O  branch.  Not 
unexpectedly,  the  Q  branch  is  not  completely  resolved,  and  appears  primarily  as  the 
broad  feature  at  123  413  cm  '.  The  S  branch  overlaps  slightly  with  the  Q  branch  and 
takes  on  a  somewhat  irregular  appearance  as  will  be  explained  later. 
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Theory  for  non-penetrating  molecular  Rydberg  states  appropriate  for  low  rota¬ 
tional  states  has  been  given  by  Herzberg  and  Jungen  (1982)  ind  by  Eyler  and  Pipkin 
(1983).  In  Hund's  coupling  case  d,  the  Hamiltonian  is  nearly  diagonal  with  the  diagonal 
elements  given  by 


EntRN  =1  -~r+R(R  + 1)  B'~aP(n, 


(1) 


The  ionisation  potential  /  is  already  given,  Ryd  for  N2  is  109  735.16  cm-1,  R  is  the 
rotational  quantum  number  of  Nj,  and  its  rotational  constant  B'  =  1.9319cm-1.  The 
polarisability  and  the  quadrupole  moment  Q  of  the  core  are  treated  here  as  indepen¬ 
dent  of  R,  as  in  the  work  of  Jungen  and  Miescher  (1969).  It  is  also  understood  that 
a  is  the  isotropic  part  of  the  polarisability,  while  the  anisotropic  part  is  essentially 
hidden  in  <?,  representing  the  effective  quadrupole  moment.  The  standard  polarisation 
formula  (Edlen  1964)  gives 


Pin.  I )  = 


3n2-I(t+l) 

2fl5(f-i)/(/  +  2l)(/+l)(/  +  i)  y 


(2) 


while  Eyler  and  Pipkin  (1983)  find 


C(l,R,  Af)  = 


3T(T-1)-4R(R  -H )/(/+!) 
2(2/  -  1)(2/  +  3)(2 R  -  1)(2R  +  3) 


(3) 


where 


Y  ~R(R  +  l)  +  Hl  +  l)-N(N  +  1).  (4) 

They  also  give  the  off-diagonal  matrix  elements  which  couple  states  of  the  same  N 
but  different  R,  however  not  odd  and  even  values  because  of  nuclear  symmetry.  These 
have  the  same  n  (and  /)  dependence  as  the  last  term  in  equation  (1).  In  diagonalising 
the  Hamiltonian,  the  value  of  Q  =  2.10  eal  found  by  Cade  (1983)  is  used  (his  value 
of  0  =  2.08  eal  is  substantially  different  from  3.12  ea%  quoted  in  Krauss  (1967)).  For 
the  7f  states,  the  shifts  of  the  eigen-energies  due  to  the  off-diagonal  matrix  elements 
is  seldom  larger  than  1cm1,  but  for  lower  values  of  n,  these  shifts  can  be  more 
substantial. 

The  initial  state  is  of  course  N2  X  l2)g  ( v  =  0,  N)  in  Hund’s  coupling  case  b,  with 
its  energy  given  by  A/ (A/  +  l)fl.  Then  the  transition  wavenumber  is  trivially  given  by 
the  energy  difference  shifted  by  aP(n,  I),  where  the  unknown  a  may  be  treated  as  a 
parameter.  However,  the  transition  strengths  are  not  readily  calculated  when  the 
lower  state  is  in  Hung’s  case  b  and  the  upper  state  is  in  case  d.  We  choose  to  transform 
the  lower  state  A-/"  into  a  superposition  of  case  d  states  preserving  parity  and 
nuclear  symmetry.  From  the  frame  transformation  theory  (Chang  and  Fano  1972), 
the  3<rg  molecular  orbital  approximated  by  the  d  component  becomes 

X?u  =  (20,  N  0|A/  -  20)4>£N .  2  +  (20,  N  0|/V  0)4>£N  +  (20,  N  0|N  +  20)<l>jJ (5) 

where  /2m2j/3mj)  is  the  Clebsch-Gordan  coefficient.  Now  the  Einstein  A 
coefficient  in  the  electric  dipole  approximation  may  be  computed  from  equation  (5) 
and  the  ni  states  4>£R  to  be 

A«N=««*(2/V'+l)(20./V()|R0rj^  ^  "2}f  (6) 


nitrogen  at  0.09  Torr  and  77  K.  Features  marked  4s  and  np  have  previously  been  identified  as  ns  o4  Flu  and  np  *nu+  respectively, 
features  marked  nf. 
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where 

A  J1  L3 

F  =  - —(21  +  1)(20, 1 0|  30)2|(3d|r  |n  f  )|2<e»  \v')2  (7) 

may  be  regarded  as  a  normalisation  factor.  The  Kronecker  delta  in  equation  (6) 
emphasises  that  the  NJ  core  remains  passive  while  the  3<rg  molecular  orbital  absorbs 
a  photon  into  the  nf  Rydberg  orbital.  However,  equation  (6)  belies  the  fact  that  an 
initial  state  is  really  a  superposition  of  three  (or  less)  core  states,  of  which  only  one 
can  contribute  to  the  A  coefficient  for  a  particular  a.  Correspondingly,  three  branches 
arise  depending  on 

f-2:  O  branch 

R  -  N  =  <  0:  O  branch 

l  +2:  S  branch. 

Since  B'  is  approximately  equal  to  B  in  N2  and  the  effect  of  the  last  term  in  equation 
(1)  is  relatively  small,  the  absorption  spectra  are  organised  into  the  above  three 
branches  with  the  Q  branch  piling  up  in  the  band  centre  and  the  O  and  S  branches 
spreading  out  with  a  spacing  of  approximately  4 B.  The  designation  of  O,  Q  and  S 
branches  to  describe  the  rotationally  resolved  photoionisation  cross  sections  in  H2 
was  advocated  by  Niehaus  and  Ruf  (1971).  Our  adoption  of  their  rotation  emphasises 
the  similarity  of  /  uncoupling  in  absorption  here  to  the  uncoupling  of  the  electron 
into  the  continuum  there. 

Further,  the  intensities  are  obtained  by  multiplying  the  A  coefficients  to  the 
population,  i.e., 

/£*  -g(2Af  +  l)exp[-Ar(iV  +  l)BAmi2N.  (8) 

In  equation  (8)  g  has  the  value  2  when  N  -  even  (nuclear  symmetry  s)  and  the  value 
1  when  S  =  odd  (a)  in  isotope  14,  which  gives  rise  to  the  alternation  of  intensities. 
At  77  K,  there  is  only  adequate  population  up  to  the  N  =  8  level.  Therefore  we  can 
expect  about  IS  lines,  which  are  easily  discernible  in  figure  2. 

In  figure  2(b),  we  show  calculated  line  positions  shifted  by  17.5  cm-1  and  the 
calculated  intensities  broadened  by  the  instrumental  line  widths.  The  lines  are  desig¬ 
nated  by  the  three  branches,  O,  Q  and  S  with  the  lower  rotational  state  given  below. 
The  subscripts  P,  Q  or  R  are  the  normal  spectroscopic  branch  notation  for  AT  -  N  =  - 1 , 
0,  or  1  respectively.  Here  they  serve  no  organising  purpose,  but  instead  cause  each 
line  to  split  into  three.  Fortunately,  in  the  O  branch,  the  0K  lines  are  substantially 
more  intense  than  the  other  two.  Consequently  they  dominate  the  spectra  from 
123  343  to  123  395  cm'1  appearing  regularly  with  a  spacing  of  about  8  cm'1  and  with 
an  alternation  of  intensities.  The  Q  branch  starts  out  at  123  405  cm'1,  but  quickly 
piles  up  at  the  broad  structure,  123  410  to  123  415  cm'1,  engulfing  the  strong  Sr(O) 
line.  Finally,  the  S  branch  assumes  a  somewhat  irregular  appearance  because  the 
intensities  here  are  more  evenly  distributed  among  the  Sp,  Sq  and  Sr  lines. 

Generally,  the  agreement  between  observed  and  calculated  line  position  is  within 
±0.5  cm'1.  Occasionally  it  is  worse  when  the  identification  is  more  doubtful.  Usually 
the  intensity  pattern  also  compares  favourably.  It  should  be  noted  that  while  the 
calculated  intensities  are  linear  the  observed  densitometer  tracings  are  not.  The  shift 
of  17.5  cm'1  in  the  calculated  spectra  implies  a  vajue  of  19  a«  for  a.  Overall,  the 
combined  position  and  intensity  matching  leave  no  doubt  of  the  correct  identification 
of  the  7f  complex. 
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An  analysis  starting  from  Hund's  case  d  has  been  used  to  extend  the  work  of  Herzberg  and  Jungen  on  the 
ig  .4/  emissions  in  H,.  A  simple  analytical  expression  for  the  line  intensities  is  presented  that  agrees  with 
their  calculations  is  about  1%.  All  of  the  experimentally  observed  lines  have  been  accounted  for  by  including 
higher  vibrational  levels  in  our  calculations. 


Recently,  a  group  of  infrared  emission  lines  in  Ht  has 
been  identified  as  the  electronic  transitions  between  two 
nonpenetrating  Rydberg  states,  namely  the  5#  and  the  4/ 
states,  by  Herzberg  and  Jungen1  (referred  to  as  HJ). 
Their  theoretical  treatment  started  from  a  Hund’s  cou¬ 
pling  case  b  basis  and  the  Hamiltonian  was  then  trans¬ 
formed  into  Hund’s  case  d.  Instead,  we  have  chosen  to 
start  from  case  d  because  the  Hamiltonian  is  already 
nearly  diagonal  in  this  basis  set,  following  the  same 
general  approach  that  was  applied  to  llu-  triplet  4tl  states 
of  Hj  by  Eyler  and  Pipkin2  (referred  to  as  EP).  After 
applying  small  corrections  towards  case  b,  the  resulting 
level  energies  are  essentially  identical  to  those  found  by 
HJ.  The  case  d  approach  is  particularly  useful  for  the 
line  intensities,  where  a  simple  analytical  expression 
gives  results  that  agree  closely  with  the  numerical  cal¬ 
culations  of  HJ. 5  This  formulation  is  nearly  identical 
to  our  treatment  of  the  analogous  states  in  atoms.  *•*  A 
number  of  the  observed  lines  were  not  previously  iden¬ 
tified.  Uy  extending  the  calculation  to  higher  vibrational 
numbers,  we  have  been  able  to  account  for  all  of  these 
lines,  indicating  that  all  of  the  observed  spectra  belong 
to  the  5*'-  4/  system. 

Iii  llund’s  coupling  case  d,  the  diagonal  term  of  the 
Hamiltonian  is  given  in  El>  as 

f.iior.lO)  =  /  -  Qjr  *  *(« ,  v)  -  aRy  (r'*  >„, 

-2Rvfy(r-3)n,4“y<r-4)Bl|(fRA'|Pf|fRfV>  .  (1) 

The  states  are  labeled  according  to  electronic  angular 
momentum  l,  core  rotation  It,  and  total  angular  momen¬ 
tum  N.  This  notation  differs  from  that  of  HJ,  but  agrees 
with  EP  and  the  commonly  accepted  usage  of  Hougen’s 
monograph.*  In  Eq.  (1)  the  first  term  is,  of  course,  the 
Rydberg  term,  where  for  Hj,  R\  =  109707.45  cm'1.  The 
next  term  is  the  rotational-vibrational  energy  of  the  H*t 
core  as  calculated  by  Hunter  el  «f.7  In  the  third  term, 
or  is  the  expectation  value  of  the  isotropic  polarizability 
of  the  core  a  =  (A/'la(A)lfti’>,  where  the  internucloar 
separation  is  denoted  by  A  to  avoid  confusion  with  the 
rotation  A.  a(A)  is  taken  from  llishop  and  Cheung*  as 
was  done  in  HJ.  The  last  term  contains  the  anisotropic 


interaction  between  the  angular  momenta  of  the  electron 
l  and  of  the  core  R.  Q  is  the  expectation  value  of  the 
quadrupole  moment  Q=(At<IQ(A)IAt;),  where  Q(A)  is 
taken  from  the  more  extensive  calculation  of  Karl  et  at.* 
rather  than  from  Bates  and  Poots10  as  was  done  in  HJ. 
The  anisotropic  polarizability  y  =  (Rv  I  or, (A)  -  orxCA)  IAt>) 
is  calculated  in  the  same  manner  as  a.  Sample  values 
of  or,  y,  and  Q  arc  3. 181  and  4.027aJ,  and  1.643en£  for 
S-l  and  a  =0,  which  compare  closely  with  those  given 
in  Table  1  of  EP. 

From  Eq.  (1),  the  energy  level  structure  of  the  5 g  (or 
4 f)  state  can  be  easily  envisioned.  All  levels  are  uni¬ 
formly  shifted  downwards  as  predicted  by  the  polariza¬ 
tion  model, 11  of  order  1  cm'1  for  the  5j/f  or  10  cm"1  for 
the  4 f.  The  main  splitting  is  due  to  the  core’s  rotational 
and  vibrational  structure  with  a  spacing  of  =>  100  cm"1, 
analogous  to  the  atomic  core  fine  structure  spacing. 
Interplay  between  A  and  /  through  the  quadrupole  inter¬ 
actions  produces  a  finer  splitting  of  the  order  =<  1  cm"1 
for  the  (or  =*  10  cm"1  for  the  4 /).  The  resultant  angu¬ 
lar  momentum  N  and,  in  turn,  the  energy  EmlKII(0)  de¬ 
pend  on  the  relative  orientation  of  l  and  R.  When  N  is 
a  maximum  (or  minimum),  /  and  R  are  parallel  (or  anti- 
paraltcl),  so  the  positively  charged  molecular  axis  lies 
in  the  plane  of  the  Rydberg  electron’s  orbit,  thereby 
lowering  the  electrostatic  energy.  In  contrast,  when  N 
is  midway,  l  and  R  are  perpendicular,  and  the  electro¬ 
static  energy  is  then  the  highest.  Of  course,  quantum 
precession  ameliorates  the  simple  classical  picture  re¬ 
sulting  in  some  out  of  the  plane  motion  and  less  electro¬ 
static  attraction  when  N  is  a  maximum.  Consequently, 
the  energy  levels  are  usually  ordered  N„„,  Nma,  Nmu  ■ 
-1,  Nmu  -  1,  etc.1  While  this  simple  picture  explains 
the  level  structure  well,  it  neglects  the  nondiagonal  in¬ 
teraction  between  different  R  levels  (subject  to  the  para/ 
ortho  restriction).  The  off-diagonal  terms  have  the 
same  form  as  Eq.  (1)  except  that  the  last  factor  must 
be  changed  to 


(IRN\  /Jj|  IR'N't  -  (-  A)J(2R<  lH2R'.  1) 

//  i  l\/R  2  R'\\N  R'  /( 

\0  0  0/\0  0  0/(2  l  R l 


(2) 


11  Present  utltlrcss:  Yale  I'nivcrsity  Physics  IM'ixirlnuMil,  New 
Haven,  C  oniieeLkcui  0051 1 . 


In  Eqs.  (1)  and  (2),  explicit  expressions  for  (r'3),,, 

(r~* )„,,  and  the  3-;  and  6 -j  symbols  were  given  in  EP. 
The  resultant  matrices  tor  each  value  of  it,  l,  »»,  and 
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TABLE  I.  Calculated  and  observed  wave  numbers  <>l  the  5f"V  lines  in  II,.  t  denotes  new 
identification  and  *  signifies  contrary  identification  to  IIJ. 


Branch  v 

®0 

Branch  „ 

“0 

2459.89 

ft, (3)  0 

2459. 50 

2500.39 

4,(2)  3 

2500.38 

2459.94 

ft,(3)  1 

2459. 50 

2502. 93 

ft, (5)  2 

2502.  89 

2460. 65 

4,(4)  2 

2459.96 

2503. 05 

ft, (2)  3 

2502.  89 

2460. 68 

ft,( 3)  2 

2459.96 

2503.63 

ft, (3)  3 

2460.69 

4,(4)  3 

2459.96 

2504. 34 

«,(1>  0 

2461.02 

4,(4)  1 

2461.44 

2504. 71 

ft, (2)  1 

2504.76 

2461.39 

4,<4t  4 

2461.44 

2504. 72 

ft, (6)  2 

2504. 76 

2461.62 

4,(4)  0 

2461.44 

2504. 90 

4,(4)  1 

2504.76 

2461.87 

ft  2  (3)  4 

2462. 14  • 

2506.  21 

R,(l)  0 

2506.36 

2462.11 

H,<3)  3 

2462. 14 

2506.52 

ft, (2)  3 

2506. 36 

2462.43 

ft,( 3)  5 

2462.47  • 

2507. 38 

4,(2)  4 

2462.48 

ft, (3)  3 

2462.47 

2508.04 
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TABLE  11.  Theoretical  line  strengths  anti  observed  intensities  according  to  the  present  assign¬ 
ments.  Brackets  tlcnulc  shared  intensities,  with  llie  square  bracket  indicating  the  strongest  line 
and  the  rouisl  bracket  weaker  lines.  0o:  blended  with  the  ulonitc  Bruekelt  alpha  line. 


Branch 

S' 

0 

1 

2 

3 

4 

5 

R0(  3) 

1.00 

12.31 

2.4 

13. 01 

0.8 

(1.41 

..  a 

A,  (4) 

3.07 

6.5 

5.5 

(2.7| 

2.3 

1.3 

0.7 

A,  (3) 

2.81 

17.01 

(7.0) 

(1.41 

(0.81 

1.3 

0.5 

R  ,(2) 

2.14 

14.71 

14.71 

1.7 

(1.6| 

see 

Ql<4) 

0.19 

(4.7) 

(4.7) 

(0.5) 

(0.5) 

(0.9) 

"* 

R2(5> 

1.44 

1.8 

(1.41 

(0.8| 

0.5 

(1.6) 

(0.9) 

R,(4) 

1.10 

1.4 

8. 

8. 

Ba 

(0.8) 

ft,  (3) 

0.82 

|1.4j 
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(0.3) 

(0.51 

A, (21 

0.0 

0.4 

(1.4) 

(0.8) 

WmWmi 

R,(11 

0.43 

10.01 

(1.4) 

(0.51 

(0.9) 

a  a  a 

*  •  * 

Aj(0) 

0.8 

0.9 

(4.7) 

(0.5| 

(0.91 

A3(5) 

3.  78 

«. 

8. 

(0.8) 

(2. 3) 

A3(4) 

2.  79 

(1.01 

10.31 

(0.  3) 

(1.0) 

(1.4) 

8„ 

A3<3> 

1.95 

10.81 

10.0) 

A3(21 

1.32 

0.5 

(2.71 

(0.6) 

0.5 

AjUI 

0.81 

(0.4| 

... 

•  .  • 

<J3(4) 

0.  9G 

[1.41 

(1.4) 

(1.41 

(1.4) 

(1.4) 

(J3(3) 

0.93 

0.5 

0.5 

(0.8) 

<J3(2I 

0.75 

0.5 

0.G 

0.8 

(2.7) 

0.7 

(-  1)K  are  diagonalized  to  yield  the  eigenenergies 
i hm  »■ 

In  conformity  with  the  HJ  notation,  the  transition  fre¬ 
quency  o  is  labeled  by  the  lower  state.  Thus, 

of/<s(lV)«]  =  £s<RW,llv  -  .  (3) 

The  intensity  is  proportional  to  the  product  of  the  line 
strength  S  and  the  Boltzmann  factor.  From  the  Einstein 
A  coefficients  given  in  EP,  the  line  strength  S  is 


S  -(2 - (—  ♦  1) /  ^  N'  R\*  .  (4) 

\N  3  1/ 

The  first  factor  is  the  statistical  weight  due  to  nuclear 
spin,  1.  e. ,  1  when  R  is  even  {para )  or  3  when  R  is  odd 
{ortho).  The  rest  of  Eq.  (4)  is  identical  to  the  analogous 
line  strengths  in  atoms.5 

Results  of  our  calculated  a  are  given  in  Table  I  along 
with  the  line  designations  and  the  observed  values  of  HJ. 


FIG.  1.  Determination  of  the  ro¬ 
tational  temperature  for  v  0 
lines.  /  given  by  Eq.  (5)  is  plot¬ 
ted  against  the  energy  of  the  he 
sublevels.  Note  the  bunching  of 
the  data  points  according  to  R 
-0,  1,  2,  and  3.  Arrows  indicate 
the  direction  of  displacement  due 
to  shared  intensities. 
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FIG.  2.  Determination  of  the  vi¬ 
brational  temperature  from  the 
ft,  (4)  lines,  7  given  by  Eq.  (5) 
is  plotted  against  the  vibrational 
energy  of  the  H$  core. 


2  4  6  8  10  12 

E  x  10 3  cm'1 

Generally,  our  calculated  values  agree  with  their  ealeu-  tunately,  the  intense  ft, (4)  lines  are  unshared  for  all  v, 

lated  ones  to  0. 1  cm"1.  In  addition,  we  have  calculated  so  a  vibrational  temperature  can  be  extracted.  Figure 

a  for  e  =  4  and  5  as  well  as  for  some  weaker  lines  with  2  shows  that  a  stra:ght  line  fits  the  data  rather  well, 

ti  =  0  to  3.  Further,  we  have  calculated  the  positions  of  where  E  is  the  vibrational  energy7  of  H*,.  From  the 

some  5/  -4d  lines.  A  few  fall  within  the  range  of  the  slope,  we  determine  a  vibrational  temperature  of  6000  K, 

HJ  spectra,  but  there  is  no  convincing  evidence  that  any  thus  continuity  the  HJ  value.  With  the  addition  of  the 

of  these  lines  have  been  observed.  It  is  possible  that  higher  vibrational  quantum  numbers,  the  task  of  iden- 

some  of  the  observed  spectra  share  intensity  with  these  tification  has  been  completed  in  that  all  experimental 
transitions.  From  Table  I,  it  is  seen  that  we  are  able  lines  are  now  labeled  and  conversely,  all  theoretically 

to  identify  all  observed  lines.  More  often  than  not,  an  strong  lines  are  matched  to  an  experimental  feature, 

experimental  feature  corresponds  to  more  than  one  cal¬ 
culated  line,  though  usually  only  one  is  expected  to  be  *n  conc^us'on>  our  extended  analysis  strengthens  the 

strong.  The  new  identifications  are  marked  t.  While  identification  of  the  spectra  between  2459  and  2534  cm" 
we  agree  with  HJ  on  the  identification  of  strong  lines,  observed  by  HJ  as  the  5# -4/  emissions  in  H,. 

we  have  made  a  few  reassignments  of  the  weaker  lines,  *n  Particular>  each  observed  line  is  assigned  to  one  or 

marked  by  asterisks.  more  ot  the  emulated  5tf-4 f  lines  from  v  =  0  to  5.  Our 

line  assignments  are  mostly  the  same  as  HJ  with  just  a 

The  task  of  line  assignment  is  greatly  facilitated  by  few  exceptions.  Thus,  there  is  no  evidence  for  the  exis- 

intensity  considerations.  In  Table  II,  we  show  the  line  fence  of  other  systems  or  resolvable  singlet-triplet 

strengths  given  by  Eq.  (4)  for  all  of  the  calculated  tran-  splittings  in  the  above  spectral  range.  However,  our 
sitions.  The  other  columns  show  the  experimental  line  extended  calculations  and  new  assignments  provide  more 
intensities  /„  from  HJ,  using  the  revised  assignments.  information  allowing  for  a  more  accurate  determination 

To  facilitate  analysis  of  the  rotational  lemperature,  the  of  the  rotational  and  vibrational  temperatures.  While 

Boltzmann  factor  was  separated  from  the  theoretical  line  file  former  is  only  moderately  well  defined  and  somewhat 

strength  from  Eq.  (4)  by  taking  the  ratio  higher  than  the  value  found  by  HJ,  the  latter  is  surpris- 

-  ingly  well  defined  and  is  identical  to  the  value  of  HJ. 

/  =  -ln(/0/S)  -E/kT  .  (5)  Finally,  this  f-uncoupling  model  including  only  long- 

In  Fig.  1,  we  plot  7  against  E  for  v  =  0.  It  fits  a  straight  range  interactions  is  shown  to  be  a  powerful  tool  for  In¬ 
line  reasonably  well,  taking  into  account  that  arrows  in-  vestigating  the  spectra  of  molecules  as  well  as  of  atoms 
dicate  the  effect  of  shared  intensity.  From  the  slope,  involving  nonpenetrating  Rydberg  states, 
we  determine  a  rotational  temperature  of  330  K  in  rough 
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TIME  RESOLVED  FTS  OF  MOLECULAR  AND  ATOMIC  INFRARED  EMISSION 
W.  Barowy  and  H.  Sakai 

Astronomy  Research  Facility,  Department  of  Physics  and  Astronomy 
University  of  Massachusetts,  Amherst,  Mass.  01003,  U.S.A. 

The  technique  of  Fourier  spectroscopy  applied  to  time  dependent 
studies  was  first  developed  by  Sakai  and  Murphy  (1,2,3).  The  technique 
is  essentially  to  measure  the  time  dependent  interferogram  signal  given 
byi 

F(x,t)  =  /B(o,t )  cos  2irox  do.  (l) 

The  spectral  signal  B(  a,  t )  is  recovered  by  applying  the  inverse 
transformation  to  the  interferogram  F(x,t)  with  the  same  ts 

B(o,t)  =  /F(x,t)  cos  2irox  dx.  (2) 

Our  scheme  is  devised  to  measure  the  time  resolved  interferogram  F(x,t) 
at  each  position  of  x  while  the  time  parameter  varies  continuously.  The 
entire  F(x,t)  are  collected  by  repeating  the  same  procedure  for  all  x 
from  sero  to  the  maximum  path  difference,  as  the  interferometer  is  held 
at  each  sample  position  and  then  stepped  to  the  next.  At  the  University 
of  Massachusetts  ,  we  are  applying  this  technique  to  the  study  of  atomic 
and  molecular  infrared  emissions. 


SOURCE 


The  source  is  a  glow  discharge  operated  in  the  central  12  meters  of 
a  30  meter  long,  1  meter  diameter  steel  tube  as  seen  in  Figure  1.  The 
discharge  is  excited  by  the  application  of  an  ac  voltage  of  up  to  1000 
volts  (peak  to  peak )  between  the  walls  of  the  tube  and  a  12  meter  long, 
cooled  aluminum  cylinder,  located  at  the  axial  center  of  the  tube.  The 
pressure  at  which  the  discharge  can  operate  reliably  is  between  50  and 
300  sullitorr.  Two  f/30  mirrors  are  situated  at  each  end  of  the  tube 
and  focused  on  each  other,  providing  an  optical  path  length  of  90 
meters.  A  KBr  lens  at  the  exit  port  of  the  tube  collimates  the  beam  in 
preparation  for  entrance  to  the  interferometer.  The  simplest  waveform 
used  for  excitation  is  a  60  Hz  sinusoid  generated  by  a  step-down 
transformer  ope rated  in  reverse.  A  pulsed  discharge  is  also  used  with 
pulses  of  both  polarities  available. 

INSTRUMENTATION 

The  detector  is  a  3  millimeter  diameter  InSb  photodiode  cooled  to 
liquid  nitrogen  temperature.  The  detector  accepts  radiation  from  the 
interferometer  at  the  focus  of  an  f/2.25  off-axis  parabolic  mirror.  The 
interferometer  is  tilt-compensated  with  corner  reflectors  and  is  scanned 
in  the  step-and-hold  mode.  Our  corner- ref lector  interferometer  is  shown 
in  Figure  2.  Sampling  of  the  time-dependent  signal  occurs  while  the 
movable  mirror  is  held  in  position  with  a  servo  mechanism  referenced  to 
the  6329A  lino  of  a  HeNe  laser.  Thus  the  unit  stepping  distance  is 
6329A. 

The  mirror  is  driven  by  a  linear  "speaker-coil"  motor  with  the 


3 

control  electronics  enslaved  to  an  LSI11/V03  computer.  The  detector 
signal  is  wide  band  amplified  (5Hz  to  10  KHz)  and  is  sampled  by  a  12  bit 
A/D  converter  at  130  microsecond  intervals,  the  cycle  time  of  the  master 
clock.  The  master  clock  is  synchronized  with  the  excitation  source  and 
the  sampling  is  "locked-in"  to  the  excitation  waveform  with  the 
phase-locked- loop  (PLL).  The  L5I11  computer  reads  in  the  sampled 
time-resolved  interferogram  F(x,t)  via  a  DRV11  parallel  1/0  interface, 
pre-processes  the  data  to  improve  the  signal-to-noise  ratio,  and  stores 
the  results  on  either  a  floppy  disk  [0X1]  or  a  magnetic  tape  [TMll]  via 
the  PDP11/20.  a  block  diagram  of  the  apparatus  is  shown  in  Figure  3. 

PREPROCESSING 

The  recording  of  the  time-resolved  interferogram  F(x,t)  is 
controlled  by  the  software  structured  as  in  the  flow  chart  shown  in 
Figure  «.  The  entire  scheme  enjoys  the  clexibility  and  power  inherent 
in  software  processing.  As  the  signal  is  repetitively  sampled,  the 
LSIll  co-adds  corresponding  points  of  each  excitation  cycle  in  real-time 
for  each  optical  path  retardation.  The  possible  number  of  coadditions 
is  virtually  unlimited;  in  reality,  however,  it  is  constrained  by  the 
holding  time  or  the  LN2  dewar.  This  mode  preserves  the  time  dependence 
of  the  signal  while  increasing  the  signal-to-noise  as  the  square  root  of 
the  number  of  coadditions.  Thus,  molecular  features  that  appear  weak  or 
out  of  phase  with  the  excitation  can  be  detected  and  the  phase  of  peak 
emission  determined  using  low  spectral  resolution.  Then  appropriate 
processing  such  as  synchronous  demodulation  can  be  programmed  for 
non-time-resolv* d,  high  resolution  studies  of  these  features. 
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However,  synchronous  demodulation  may  distort  a  spectrum  in  a 
modulating  source  due  to  different  excitation  and  de-excitation  rates 
and  phases  of  individual  spectral  features  (4).  In  this  case,  where 
interest  is  in  the  entire  spectrum,  an  asynchronous  sum  is  done,  which 
is  simply  the  integral  of  the  modulus  of  the  A.C.  component  of  the 
waveform.  In  this  manner  all  frequency  components  of  the  signal  are 
summed,  independent  of  phase,  and  the  contributions  of  each  spectral 
element  are  maximized.  60  Hz  pick-up,  which  is  independent  of  path 
difference,  adds  a  D.C.  offset  to  the  interferogram  in  this  type  of 
average.  The  offset  can  be  calculated  however,  by  sampling  the  signal 
with  a  shunt  resistance  placed  across  the  detector  output.  This  can 
then  be  subtracted  later  from  the  interferogram.  At  present, 
preprocessing  of  non-time-resolved  spectra  allow  integration  times  that 
are  constrained  solely  by  the  dewar  holding  time. 

Instabilities,  such  as  arcs  in  the  discharge,  are  compensated  by 
intelligent  data  acquisition  routines.  This  is  especially  important  in 
Fourier  Spectroscopy,  where  an  arc  could  cause  a  spike  in  the 
interferogram,  thereby  distorting  the  entire  spectrum. 

By  a  short  study  of  arcing  phenomena,  we  have  found  that  arcs  tend 
to  appear  at  specific  times  in  the  excitation  cycle.  When  an  arc 
occurs,  it  is  apparent  in  both  the  waveform  of  the  disharge  current  as 
well  as  the  detector  signal.  Implementation  of  "smart"  routines  with 
time  resolution  allow  the  computer  to  identify  these  anomalies  within 
certain  toleran>  «*s,  and  to  reject  them  if  they  occur. 


POSTPROCESSING 


The  spectral  recovery  process  is  performed  as  indicated  in  Figure 
4.  The  first  step  in  rearranging  F(x,t)  requires  a  very  large  working 
space  in  order  to  perform  it  with  high  efficiency.  For  example,  a  study 
with  spectral  resolution  of  1  wavenumber  requires  over  one  million  data 
to  be  shuffled.  Currently  we  use  our  University  central  site  CDC  Cyber 
system  for  the  rearranging  process  as  well  as  the  rest  of  the  spectral 
recovery . 

PRELIMINARY  RESULTS 

Some  typical  data  for  nitrogen  and  oxygen  in  the  60  Hz  discharge 
are  shown  in  Figures  S  and  6,  respectively.  The  time  evolution  flows 
from  the  top  to  the  bottom  of  the  figures.  For  clarity,  the  full  time 
resolution  is  not  displayed,  rather,  the  presented  spectral  plots  occur 
approximately  2  milliseconds  apart.  The  spectral  resolution  is  4  wave 
numbers.  The  first  plot  of  each  series  corresponds  to  a  point  in  time 
when  the  voltage  of  the  central  electrode  is  positive  with  respect  to 
the  walls  of  the  tube  and  near  its  peak  value. 

Several  features  in  the  development  of  the  spectrum  of  nitrogen 
that  are  immediately  apparent  are  the  B3n-A3E  transitions  in  Figure 
5(e-h)  between  0000  and  7500  wave  numbers.  Atomic  lines  of  nitrogen  and 
oxygen  in  Figure  5(e-g)  and  Figure  6(e-g)  show  preferentially  during 
only  specific  stages  of  the  discharge  cycle. 

We  believe  that  the  cylindrical  geometry  of  the  source  leads  to  a 
significant  di t  terence  in  the  distribution  of  electron  energies  with 
respect  to  the  polarity  of  the  applied  voltage.  In  general,  the  cross 
sections  for  o.ectron  impact  excitation  from  the  ground  state  of  the 
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molecule  to  an  upper  state  depend  on  the  identity  of  the  excited  state 
as  well  as  the  incident  electron  energy  (S).  As  a  result,  individual 
electronic  states  of  a  molecule  manifest  themselves  in  the  discharge  at 
different  stages  as  seen  in  Figures  5  and  6.  This  is  complicated  by  the 
extreme  variety  in  lifetimes  of  these  electronic  states  (6,7)  and  by  the 
role  of  collisions,  which  lead  to  such  phenomena  as  intra-system 
cascading . 

The  process  of  unraveling  these  effects  will  be  facilitated  with 
the  pulsed  discharge,  where  the  external  excitation  can  be  applied  or 
removed  rapidly,  and  the  evolution  of  excited  states  followed  by 
examining  the  infrared  spectrum.  At  our  present  operating  pressures, 
many  of  the  infrared  transitions  of  atmospheric  molecules  have  radiative 
decay  rates  comparable  to  the  mean  inter-molecular  collision  rate, 
thereby  providing  a  good  probe  for  exploring  the  excitation  and 
relaxation  phenomena. 
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Figure 

1  Schematic  representation  of  the  discharge  source. 

2  Corner-Reflector  Interferometer 

3  Block  Diagram  of  the  Apparatus 

U  Flow  Chart  for  Data  Processing 

5  The  time  evolution  of  the  spectrum  of  molecular 

nitrogen  in  the  modulated  glow  discharge. 
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The  time  evolution  of  the  spectrum  of  molecular 
oxygen,  showing  several  atomic  lines  appearing 
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The  infrared  radiative  process1  through  the  lower  part  of 
atmosphere  is  a  composite  effect  of  the  local  characteristics,  which  are 
principally  controlled  by  three  parameters,  the  local  abundance  of  the 
molecules  responsible  for  the  radiative  process,  the  local  temperature 
and  the  local  pressure.  The  molecules  in  the  lower  layers  of  the 
atmosphere  are  generally  in  thermal  equilibrium  with  the  local 
surrounding.  The  third  parameter,  the  local  pressure,  plays  its 
effective  role  only  in  the  atmosphere  below  the  tropopause,  while  the 
other  two  in  their  role  for  controlling  the  infrared  radiative  process 
are  important  throughout  the  entire  low  atmosphere  including  the 
stratosphere .  Under  the  local  thermal  equilibrium  the  radiative 
transfer  characteristics  in  infrared  are  determined  by  integrating  the 
local  characteristics  in  a  straightforward  computation,  once  these 
parameters  sure  known  along  the  path.2  Thus  throughout  the  atmosphere  up 
to  the  stratosphere  ,  these  parameters  are  important  in  their  critical 
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role  of  controlling  the  overall  budget  of  the  incoming  and  the  outgoing 


thermal  radiation.  At  present,  our  Knowledge  of  these  parameters  in  the 
stratosphere  is  extremely  scarce , as  a  result  of  limited  observed  data 
existing  at  present.  In  hoping  to  improve  our  Knowledge  of  these 
parameters  in  the  stratosphere,  we  designed  an  experiment  to  collect  the 
infrared  radiative  transfer  characteristics  in  the  stratosphere  by 
observing  the  infrared  photons  emitted  along  the  horizontal  line  of 

sight.  A  cryogenically-cooled  Fourier  interferometer  spectrometer 

3-5 

mounted  on  a  balloon-borne  platform  was  used  for  the  experiment .  Tne 
balloon  flight  tooK  place  on  October  7,  1981,  at  Holloman  AFB,  New 
Mexico,  and  obtained  the  spectral  data  taKen  at  the  balloon  altitude  of 
28500  m  for  more  than  two  hours. 

In  the  past  ,  the  mainstream  of  the  stratospheric  spectral 
study  carried  out  by  using  the  balloon-born  platform  was  the  absorption 
scheme,  i.e.,  by  means  of  the  solar  spectra.6  A  problem  central  to  the 
spectral  study  of  the  stratosphere  ,  a  long  absorption  path  necessary 
for  the  detection  of  various  atmospheric  species,  is  not  solved  in  a 
satisfactory  degree  with  the  solar-spectra  scheme,  since  the  absorption 
path  is  locKed  onto  the  solar  position,  its  extreme  length  is  available 
at  the  time  of  sunrise  or  sunset.  The  time  available  for  such 
measurement  is  limited  to  a  short  time.  In  addition  to  this  problem  , 
the  scheme  does  not  allow  a  diurnal  study  of  some  molecules  which  are 
suspected  of  concentration  variation  with  respect  to  the  solar  position. 
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The  spectrometric  study  on  the  atmospheric  emission  was 
conducted  using  the  satellite-borne  spectrometer.7  The  measurements  were 
to  observe  a  cumulative  effect  of  the  atmosphere  along  the  vertical  line 
of  sight.  The  individual  contribution  of  a  thin  atmospheric  layer  to 
the  overall  radiative  transfer  characteristics  was  determined  by  an 
elaborate  analysis  scheme,  which  provided  no  direct  grip  on  the  local 
characteristics  of  the  radiative  trams fer . 

Our  study  is  intended  to  supplement  these  two  spectrometric 
studies  by  measuring  the  atmospheric  emission  spectra  along  the 
horizontal  line  of  sight.  The  measurements  are  free  from  the 
limitations  imposed  on  the  solar  spectrum  scheme.  They  would  be  done  at 
any  time  of  day,  any  season,  and  at  any  geographical  location.  In 
addition  ,  the  measurement  scheme  is  conceptually  simple  because  no 
elaborate  solar  tracking  system  is  required.  The  data  collection 
extends  from  the  beginning  to  the  end  of  the  balloon  flight.  Our  scheme 
obtains  a  long  atmospheric  path  in  a  thin  slice  of  the  atmosphere  along 
the  horizontal  line  of  sight. 

The  radiance  level  expected  for  the  atmospheric  emission  is 
characterized  by  the  blackbody  radiance  of  the  stratospheric  temperature 
at  approximately  220  K.  The  peak  radiance  level  corresponds  to  an 
optically  thick  layer  of  the  emitting  molecules  .  in  general  the 
radiance  level  is  expected  less  than  the  blackbody  radiance  at  220  K. 
Even  though  the  long  path  along  the  horizontal  line  of  sight  would 
optimize  the  obtained  infrared  radiance  level,  it  would  remain  extremely 
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low,  thereby  requiring  a  superb  sensitivity  to  the  spectrometer.  The 

instrumentation  problem  imposes  a  challenge  for  achieving  the  ultimate 

8-10 

sensitivity  of  Fourier  spectroscopy,  since  the  incoming  photon  density 
is  lower  than  the  blackbody  radiance  density  of  220  K,  the 
interferometer  spectrometer  itself  is  required  to  operate  at  a 
temperature  below  the  stratospheric  temperature  of  220  K. 

Interferometer 

Spectroscopy  involved  in  this  experiment  requires  a  high 
sensitivity,  a  high  spectral  resolution,  and  a  wide  spectral  coverage. 
The  technique  of  Fourier  spectroscopy  is  only  the  scheme  which  satisfies 
these  requirements.  The  interferometer  which  is  a  central  part  of  the 
instrumentation  must  function  adequately  at  a  temperature  much  below  the 
stratospheric  temperature.  With  consideration  of  various  factors 

involved  in  operation  of  the  instrument,  we  decided  to  operate  the 
interferometer  at  77  K.  The  mechanical  stability  over  a  wide 

temperature  range  is  a  paramount  importance  required  to  the 
interferometer.  James  Pritchard  took  chief  responsibility  for  its 
design  and  fabrication.  All  components  of  the  interferometer  were 
fabricated  with  A 2  steel.  They  underwent  numerous  annealing  cycles  to 
reduce  their  thermal  stress  until  they  showed  a  predictable  thermal 
behavior  in  these  annealing  cycles  between  300  K  and  77  K.  The  entire 

interferometer  was  assembled  on  a  solid  base  plate,  as  shown  in  Figure 

l.  We  adapted  the  cat's  eye  retroreflecting  optical  system  for  the 
compensation  of  the  wavefront  tilt  so  as  to  make  the  two  interfering 
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beans  always  stay  in  alignment.  The  ways  which  support  the  movable 
carriage  provided  a  smooth  motion  to  it  at  77  K.  The  ball-bearing  slide 
together  with  dry  lubricant  was  found  to  satisfy  the  mechanical 
requirement  even  at  this  severe  temperature  condition.  The  movable 
mirror  was  driven  by  a  speaker-coil  type  magnetic  motor  and  its  position 
was  monitored  by  both  a  mechanical  trams  dues  r  and  the  interference 
fringes  formed  by  the  HeNe  cw  lamer  beaun.  The  beamsplitter  of  a  Ge  film 
coated  on  a  KC1  substrate  warn  paired  with  a  compensator  plate  of  a 
similar  optical,  path,  and  they  were  mounted  jointly  to  the  supporting 
frame,  using  a  spring-loaded  retaining  ring  which  served  as  a  cushion  to 
absorb  a  thermal  stress  associated  with  the  temperature  cycle.  The 
essential  part  of  our  instrument  ,  the  interferometer  spectrometer 
together  with  the  detector  assembly,  was  housed  in  a  cryogenic  chamber 
which  maintained  its  interior  temperature  at  77  K.  The  detector  warn  a 
copper-doped  germanium  photoconduct ive  element  operating  at  the  liquid 
helium  temperature.  The  detector  and  its  dewar  were,  thus,  encased  in 
the  77  K  chamber.  Figure  2  shows  the  interior  of  this  chaunber .  The 
mounting  plate  seen  at  the  middle  of  the  chamber  wa is  supported  by  four 
fiberglass-reinforced  epoxy  pillaurs .  The  interferometer  and  the 
detector  assembly  mounted  on  the  base  plate  were  secured  onto  this 
mounting  plate.  The  liquid  nitrogen  reservoir  seen  in  the  lower  section 
of  the  chamber  held  enough  liquid  nitrogen  to  maintain  the  interior 
temperature  at  77  K  for  10  hours  or  longer. 


The  interferometer  scanning  was  monitored  by  the  BeNe  laser 
line  at  6329A.  The  laser  beam  was  introduced  into  the  chamber  interior 
through  an  optical,  fiber  from  the  laser  unit  mounted  outside  the 
chamber.  This  arrangement  worked  almost  adequately,  even  though  with 
some  mystery,  during  some  time  period  in  the  balloon  flight  the  laser 
fringe  signal  modulation  diminished  to  a  low  figure.  Fortunately,  the 
modulation  returned  to  a  healthy  value  during  the  last  half  of  the 
observation  period. 

The  interferometer  drive  system  was  designed  to  reproduce  its 
scanning  motion  with  a  high  degree  of  accuracy.  Its  starting  position 
was  to  occur  within  l.Oum  at  a  distance  of  250wm  from  the  central 
maximum  position  of  the  interferogram.  The  scanning  was  to  change  its 
direction  approximately  at  a  distance  of  10  cm  on  the  other  side  of  the 
central  maximum  position.  The  mechanical  accuracy  of  this  degree  was 
required  for  the  interferometer  scanning  because  of  two  maun  reasons: 
( 1 )  no  optical  means  of  detecting  the  central  maximum  position  could  be 
implemented  in  our  setup  because  they  would  cause  a  devastating  effect 
on  the  spectral  sensitivity  of  the  instrument;  and  (2)  the  detector 
electronics  scheme  was  designed  to  make  use  of  the  reproducible 
interferometer  scanning.  One  complete  cycle  of  the  interferometer 
scanning  consisted  of  20-seconds  forward  motion  during  the  interferogram 
measurement  phase  and  10-seconds  reverse  motion.  The  detector  amplifier 
was  programmed  to  switch  its  gain  to  a  higher  setting  at  the  800th  data 
point  from  the  starting  position.  Any  malfunction  in  the  scanning  would 
cause  an  irreparable  damage  to  the  interferogram  data  collection,  since 
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the  central  maximum  and  the  consequent  gain  switching  would  occur  at  a 
wrong  place.  We  were  able  to  maintain  the  interferometer  scan 
satisfactorily  during  the  entire  balloon  flight. 


we  can  write  the  signal  intensity  s  reaching  the  detector  by 


S  =  BAft(Ao)Tn 


using  the  parameters  defined  belowt 

B  -  signal  radiance  per  a  unit  cross-sectional  area  of 
the  incident  beam,  a  unit  solid  angle,  a  unit  wavenumber 
interval,  and  a  unit  time; 

A a  ”  spectral  resolution; 

A  -  aperture  area  of  the  interferometer  mirror; 
ft  -  solid  angle  subtended  by  the  detector; 

T  -  time  period  covered  in  the  measurement;  and 
n  “  optical  transmission  efficiency. 


The  solid  angle  ft  subtended  by  a  detector  circularly  shaped 
is  determined  by  9 


ft  =  2tt/R, 


where  R  is  the  resolving  power  at  the  highest  wavenumber  end  of  the 
observed  spectral  bandpass.  The  amount  of  noise  generated  by  the 
detector  during  the  same  time  interval  is  given  by 


/Ad  /r 
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where  D*  is  the  normalized  detectivity  of  the  detector,  and  Ad  is  the 


detector’s  area. 


The  measurement  sensitivity  cam  be  conveniently  specified  by 
the  noise  equivalent  spectral  radiance 


NESR 


/r 

D*AQn(Ao) 


which  is  equivalent  to  the  spectral  radiance  giving  the  S/N  of  unity. 
The  detector’s  noise  is  controlled  by  two  factors,  the  thermal  noise  of 
the  detector  at  its  operating  temperature  and  the  fluctuation  of  the 
background  photon  stream.  With  reducing  the  spectrometer ' s  operating 
temperature,  the  latter  contribution  to  the  detector  noise  should  be 
substantially  reduced.  The  detector  sensitivity  is,  therefore,  expected 
to  show  an  improvement . 


Electronics 


The  inter ferogram  signal  level  was  expected  to  vary  a  large 
amount  by  more  than  one  order  of  a  magnitude  between  the  maxi mum  at  the 
ground  level  and  the  minimum  at  the  ceiling  altitude  of  approximately  30 
km.  The  detector  amplifier  was  designed  to  accommodate  this  large 
signal  variation  by  implementing  an  automatic  gain  adjustment  scheme. 
The  lower  gain  setting  of  the  amplifier  in  the  central  m»ir-inmm  region 
was  adjusted  by  degree  of  the  previous  central  maximum  modulation.  As 
mentioned  earlier,  the  amplifier  was  programmed  to  switch  its  gain 
setting  to  a  higher  value  by  the  ratio  of  5.64  at  the  800th  data  point 
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from  the  starting  position  of  the  interferogram  scanning. 

The  digitized  interferogram  data  and  other  pertinent  signals 
for  the  experiment  were  transmitted  to  the  ground  station  through  a 
telemetry  radio  link.  A  schematic  block  diagram  shown  in  Fig.  3 
illustrates  the  path  of  all  electronic  signals  included  in  our 
experiment.  A  crucial  requirement  for  the  data  transmission  was  to 
satisfy  the  condition  that  in  Fourier  spectroscopy  every  interferogram 
data  point,  equally  critical  to  the  spectral  recovery,  should  not  be 
lost  in  any  circumstances .  The  data  electronics  were  designed  to 
satisfy  this  requirement.  The  interferometer  scanning  was  set  at  a  rate 
of  4500  interferogram  data  per  second  for  20  seconds  in  a  single  scan. 
The  data  were  made  into  a  16-bit  word  format,  a  12-bit  interferogram 
data,  and  an  additional  4-bit  interferogram  status.  The  telemetry 
signal  was  made  into  a  total  of  72-bit  string,  which  consisted  of  a 
20-bit  synchronization  code,  a  16-bit  interferogram  data,  a  4-bit  frame 
count,  and  four  other  data.  It  was  transmitted  at  book  bits  per  second, 
or  at  approximately  11. IK  telemetry  frames  per  second.  Thus,  the 
telemetry  data  were  transmitted  at  a  rate  approximately  2.5  times  faster 
than  the  interferogram  data  acquisition,  thereby  accommodating  the 
non-uniform  interferometer  scanning  speed. 


Flight  Package 


The  flight  package  shown  in  Fig.  4  was  assembled  by  the 
University  of  Denver.  The  infrared  atmospheric  emission  signal  was 
admitted  to  the  cold  chamber  interior  via  a  ZnSe  window  which  remained 
at  ambient  temperature.  An  assembly  which  consists  of  a  plane  mirror 
intercepting  the  interferometer's  line  of  sight  to  direct  toward  a 
vertical  direction  and  a  blackbody  source  was  implemented  for  the 
purpose  of  providing  the  down- looking  observation  and  the  on-board 
radiance  level  calibration.  This  assembly  was  mounted  just  outside  the 
ZnSe  window  with  a  trolley  structure  which  could  place  the  assembly  at 
three  positions  by  the  command  sent  from  the  ground  station.  The  first 
position  called  for  a  clear  view  to  allow  interferometer's  line  of  sight 
at  the  horizontal  direction.  Nothing  obstructed  its  view.  The  second 
was  to  place  the  blackbody  source  in  the  field  of  view  for  accommodating 
the  radiance  level  calibration,  and  the  third  was  to  direct  the  line  of 
sight  toward  the  downward  direction  by  placing  the  intercepting  mirror. 
In  addition  to  this  mirror-blackbody  arrangement ,  the  entire  flight 
package  was  set  to  tilt  by  the  command  given  from  the  ground,  with  this 
arrangement  the  interferometer's  line  of  sight  could  be  raised  up  to  10 
degrees  from  the  horizon. 

Balloon  Flight 

The  flight  profile,  given  in  terms  of  the  balloon  altitude  as 
a  function  of  time,  was  monitored  using  radar  on  the  ground.  It  is 
shown  in  Fig.  5  with  the  local  temperature  measured  with  a  thermistor 
bolometer  hung  below  the  gondola.  The  measured  value  of  the  total 


radiance  level  is  also  shown  in  the  figure.  Even  though  the  laser 
signal  failure  which  occurred  during  some  part  of  the  flight  resulted  in 
some  loss  of  the  data,  the  interferometer  produced  the  analyzable  data 
for  more  than  two  hours  at  the  maximum  balloon  altitude.  The  spectral 
data  obtained  during  the  flight  consisted  of  three  parts t  (l)  the  data 
along  the  horizontal  line  of  sight;  (2)  the  data  along  the  vertical  line 
of  sight;  and  (3)  the  radiance  calibration  spectra  obtained  using  an 
on-board  blackbody  source  which  was  presumably  in  thermal  equilibrium 
with  the  ambient .  As  seen  in  the  thermistor  bolometer  data,  the  ambient 
temperature  at  the  altitude  of  28000  m  'v  28SOO  m  fluctuated  between  270 
K  and  250K,  closely  correlating  with  the  balloon  motion.  The  data 
indicated  that  the  stratospheric  ambient  temperature  measured  with  the 
on-board  thermistor  was  considerably  higher  than  that  of  the  model 
atmosphere,  and  that  its  fluctuation  was  unexpectedly  large. 

Our  balloon  flight  was  accompanied  by  two  radio-sonde  flights. 
The  temperature  measured  by  the  radio-sonde  sensor  was  at  most  230K,  as 
seen  in  Fig.  6,  at  this  altitude  range,  much  lower  than  the  temperature 
measured  with  the  on-board  thermistor.  The  on-board  instrument  consumed 
a  power  of  approximately  250W  ^  30ow.  Since  there  was  no  energy  sink 
provided  in  the  flight  package,  the  consumed  energy  by  the  on-board 
instruments  was  eventually  removed  to  the  ambient  atmosphere  through 
thermal  radiation.  Thus,  the.  temperature  of  the  instrument  package  had 
to  remain  considerably  higher  than  the  local  surrounding.  A  difference 
of  20  'v*  3 OK  between  the  temperature  registered  by  our  balloon-borne 
thermistor  and  the  radio-sonde  sensor  w as  found  within  the  estimate 


given  based  on  the  general  configuration  of  the  flight  instrument . 
Calibration  of  the  Radiance  Level 

Two  sets  of  spectra  were  taken  using  the  on-board  blackbody 
source ,  one  at  17 : 30  GMT  and  another  between  18 1 45  and  19  s  30  GMT .  Both 
sets  were  taken  at  the  balloon  ceiling  altitude  of  28000  m.  Their 
spectral  characteristics  as  well  as  their  radiance  levels  were  found 
very  similar,  we  estimated  these  spectra  to  be  those  of  a  255K 
blackbody  radiance,  and  proceeded  to  calibrate  all  the  measured  spectral 
data,  horizontal-looking  and  down- looking,  accordingly.  The  calibrated 
down- looking  spectra  showed  the  radiance  level  at  the  center  of  the  C02 
15  y  band  corresponding  to  215K  'v  220K  and  that  of  the  window  region  to 
the  ground  temperature  of  285K.  The  radiance  level  obtained  for  the 
down-looking  spectra  was,  therefore,  consistent  with  the  atmospheric 
radiative  transfer  expected  for  this  spectral  region. 

The  radiance  level  of  the  interferogram  signal  was  estimated 
with  respect  to  the  noise  level  contained  in  it.  Since  no  information 
on  the  gain  factor  of  the  detector  amplifier  was  transmitted  to  the 
ground,  it  had  to  be  determined  from  the  recorded  interferogram  data. 
The  noise  level  in  the  interferogram  data  increased  or  decreased 
whenever  the  detector  amplifier  switched  its  gain  setting.  Since  we  had 
no  reason  to  believe  that  the  absolute  level  of  noise  in  the 
interferogram  signal  varied  from  time  to  time,  we  determined  the  gain  of 
the  detector  amplifier  using  the  noise  level  observed  in  the  recorded 


interferogram  signal.  We  observed  the  noise  level  stepped  in  accordance 
with  the  step  in  the  amplifier  gain  setting.  Therefore,  our  procedure 
for  estimating  the  detector  asplifier  gain  was  self-consistent.  The 


total  radiance  level  received  by  the  spectrometer  was  proportional  to 
the  central  maximum  modulation  and  it  was  estimated  by  multiplying  the 
central  maximum  modulation  with  the  gain  factor  which  was  determined 
from  the  noise  level  observed  in  the  interferogram  very  far  from  the 
central  maximum. 

An  exception  to  the  procedure  for  the  radiance  calibration  had 
to  be  observed  for  the  data  taken  after  19:30  GMT,  when  the  blackbody 
source  was  removed  from  the  field  of  view  at  the  last  time  in  the 
flight.  The  radiance  level  determined  using  the  procedure  above  showed 
a  dramatic  increase  for  the  horizontal-looking  spectra  taken  after  19:30 
GMT,  while  the  balloon  altitude  and  the  measured  ambient  temperature 
remained  unchanged  as  before.  At  that  time,  the  liquid  helium  level  in 
the  detector  dewar  became  definitely  low  and  the  detector  temperature 
started  to  rise.  We  reasoned  that  the  apparent  increase  of  the  radiance 
level  after  19:30  GMT,  shown  by  a  solid  curve,  was  resulted  by  an 
improved  noise  figure  of  the  detector,  rather  than  by  an  increase  in  the 
incident  radiance  level.  The  radiance  level  probably  remained  at  the 
same  level,  as  indicated  by  a  dotted  curve.  Anyway,  the  data  obtained 
after  19:30  GMT  certainly  had  a  larger  degree  of  ambiguity  than  the  rest 
because  the  radiance  level  calibration  was  not  repeated  after  the 
observation  of  these  spectra.  Our  discussion  on  the  obtained  data  will 
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be  made  excluding  the  data  collected  after  19:30  GMT. 


The  spectral  response  of  our  spectrometer  was  determined  using 
the  blackbody  calibration  spectrum.  As  described  above,  the  calibration 
source  was  assumed  at  255  K.  The  obtained  blackbody  spectrum  and  the 
determined  radiance  level  given  as  a  function  of  wavenumber  are  shown  in 
Figure  7.  The  obtained  curves  were  used  to  determine  a  relative  value 
of  the  radiance  level  at  each  wavenumber  position.  The  absolute  value 
of  the  radiance  level  was  then  determined  by  normalizing  the  integrated 
relative  radiance  value  to  the  total  radiance  level  obtained  using  the 
procedure  described  above. 

Recovery  of  the  Spectral  Data 

The  telemetry  signals  recorded  on  magnetic  tapes  were  brought 
back  to  the  University  of  Massachusetts  for  retrieval  of  the  data 
obtained  during  the  flight.  First  the  telemetry  tapes  were  played  back 
and  the  inter ferograa  data  were  extracted  using  the  PDP11  computer  with 
a  PCM  telemetry  signal  decoder.  Then  the  recovered  raw  interferogram 
data  were  recorded  on  the  7-track  digital  magnetic  tapes.  The  rest  of 
the  numerical  process  was  done  on  the  CDC  Cyber  computer  system  at  the 
University  of  Massachusetts.  Since  some  of  the  data  were  sampled  with 
faulty  sampling  intervals  caused  by  insufficient  laser  interference 
modulation,  the  interferogram  data  were  inspected  one  by  one  prior  to 
application  of  the  numerical  Fourier  transformation  process.  The 
interactive  feature  developed  by  us  on  the  CDC  system  provided  an 
extremely  efficient  inspection  on  the  raw  data  displayed  on  the  CRT 
terminal.  The  phase  correction  operation  which  followed  the 


interferogram  inspection  was  applied  individually  to  the  examined 
inter ferogram  data.  The  spectral  recovery  process  was  carried  out  with 
the  sine  instrument  function  of  a  basic  sampling  interval  of  0.24  cm**1. 
Even  though  the  interferogram  data  extended  more  than  the  optical  path 
difference  of  4.2  cm,  a  limited  signal-to-noise  ratio  in  the  data 
necessitated  a  premature  termination  of  the  data  processing  at 
approximately  2.07  cm  of  the  optical  path  difference  [corresponding  to 
the  resolution  figure  of  0.48  cm-1  ]. 


The  Obtained  Spectral  Data 

The  major  atmospheric  bands  observable  in  our  measurement  are 
co  2  (  Av2  -  1 )  transitions } 1  All  bands  contain  a  strong  Q  branch .  Table 
I  lists  the  band  center  for  the  observable  transitions.  The  Q  branch  of 
the  (01101-00001)  band  of  l2C1602  isotope  should  show  a  strong 
saturation  even  at  our  balloon  altitude  .  The  radiance  level  at  the 
center  of  the  Q  branch  serves  as  a  temperature  indicator  of  the  air 
surrounding  the  balloon  .  With  the  assumption  that  the  C02 
concentration  in  the  stratosphere  is  given  by  the  mil-established 
mixing  ratio  of  3.3  x  10 “4  ,  our  observation  along  the  horizontal  line 
of  sight  was  expected  to  be  made  with  the  C02  column  density  of 
approximately  5.0  x  lo2iilolecules/cm2  . 


The  emission  spectra  observed  along  the  vertical  line  of  sight 
is  presented  in  Figure  8.  All  data  observed  along  the  vertical  line  of 
sight  differed  very  little  among  them  in  their  spectral  characteristics 
as  well  as  of  their  radiance  level.  The  source  radiance  level  in  the 
675  cm  region  was  found  to  be  approximately  4.0  x  10” 6  W/cm2  str  cm”1, 

the  blackbody  radiance  level  of  210^220  K.  The  level  in  the  800  cm”1 
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was  1.1  x  10  W/cm  str  cm  ,  consistent  with  the  concurrently  measured 

ground  temperature  of  285  K.  The  theoretical  spectrum  calculated  using 

12  11 

the  FASC0D1B  code  with  the  AFGL  atmospheric  line  parameters  is  shown  in 

Figure  9.  No  significant  differences  were  observed  between  the  measured 
and  calculated  data. 

The  data  observed  along  the  horizontal  line  of  sight  showed 
their  spectral  characteristics  very  little  different  from  each  other. 
The  difference  observed  in  their  radiance  level  is  clearly  indicated  in 
the  data  of  Figures  10  and  11,  which  were  taken  at  18:15  GMT  and  18:45 

GMT.  A  continuum  feature  which  amounted  to  the  radiance  level  of  2  a.  3x 

•  6  2 

10  w/cm  str  cm  at  the  600  cm  region,  is  probably  due  to  the 
emission  of  the  ZnSe  window.  A  large  fluctuation  in  the  radiance  level 
at  the  675  cm-1  region  indicates  a  rather  large  temperature  fluctuation 
in  the  ambient  atmosphere. 

A  theoretical  spectrum  shown  in  Figure  12  was  generated  using 
the  FASC0D1B  with  the  C02  column  density  of  1.2xl021  molecules/cm2  at 
220  K.  The  theoretical  and  the  observed  disagree  by  a  large  amount . 
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Interferometer  assembled  on  a  base  plate. 
Interior  of  Cryogenic  Chamber. 

Front  View  of  the  Entire  Package. 
Schematic  Diagram  of  Instrumentation. 
Profile  of  Flight,  October  1981. 
Radio-Sonde  Records. 


Blackbody  Source  Calibration  Curve. 

Spectrum  observed  along  the  vertical  line  of  sight. 

Synthetic  spectrum  computed  for  the  vertical  line  of  sight. 
Spectrum  observed  along  the  horizontal  line  of  sight:  18:15  GMT. 
Spectrum  observed  along  the  horizontal  line  of  sight:  l8:U5  GMT. 

Synthetic  spectrum  computed  for  the  horizontal  line  of  sight: 

21  2 
CO^  column  density  =  1.2  x  10  molecules/cm  , 


temperature  =  220°K. 
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